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Abstract 

Feynman's path amplitude formulation of quantum mechanics is used to analyse 
the production of charged leptons from charged current weak interaction processes. 
For neutrino induced reactions the interference effects predicted are usually called 
'neutrino oscillations'. Similar effects in the detection of muons from pion decay are 
here termed 'muon oscillations'. Processes considered include pion decay (at rest 
and in flight), and muon decay and nuclear /?-decay at rest. In all cases studied, 
a neutrino oscillation phase different from the conventionally used one is found. A 
concise critical review is made of previous treatments of the quantum mechanics of 
neutrino and muon oscillations. 
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1 Introduction 



The quantum mechanical description of neutrino oscillations [1, 2] has been the sub- 
ject of much discussion and debate in the recent literature. The 'standard' oscillation 
formula [3], yielding an oscillation phase 1 , at distance L from the neutrino source, be- 
tween neutrinos, of mass m\ and m 2 and momentum P, of: 2 

= M (i.i) 

is derived on the assumption of equal momentum and equal production times of the two 
neutrino mass eigenstates. Other authors have proposed, instead, equal energies [4] or 
velocities [5] at production, confirming, in both cases, the result of the standard for- 
mula. The latter reference claims, however, that the standard expression for 0i2 should 
be multiplied by a factor of two in the case of the equal energy or equal momentum 
hypotheses when different production times are allowed for the two mass eigenstates. 
However, the equal momentum, energy or velocity assumptions are all incompatible with 
energy-momentum conservation in the neutrino production process [6]. In two recent 
calculations [7, 8] a covariant formalism was used in which exact energy-momentum con- 
servation was imposed. These calculations used the invariant Feynman propagator [9] to 
describe the space-time evolution of the neutrino mass eigenstates. In Ref. [8] a formula 
for the neutrino oscillation phase differing by a factor of 9.9 from Eqn(l.l) was found for 
the case of pion decay at rest, and it was predicted that correlated spatial oscillations in 
the detection probability of neutrinos and the recoiling decay muons could be observed. 
However, the author of Ref. [7] as well as others [10, 11] claimed that muon oscillations 
would either be completely suppressed, or essentially impossible to observe. 

The present paper calculates the probabilities of oscillation of neutrinos and muons 
produced by pions decaying both at rest and in flight, as well as the probabilities of neu- 
trino oscillation following muon decay or (3— decay of a nucleus at rest. The calculations, 
which are fully covariant, are based on Feynman's reformulation of quantum mechan- 
ics [12] in terms of interfering amplitudes associated with classical space-time particle 
trajectories. The essential interpretational formula of this approach 3 , though motivated 
by the seminal paper of Dirac on the Lagrangian formulation of quantum mechanics [13], 
and much developed later in the work of Feynman and other authors [14], was actually 
already given by Heisenberg in 1930 4 [15]. The application of the path amplitude formal- 
ism to neutrino or muon oscillations is particularly staight forward, since, in the covariant 
formulation of quantum mechanics, energy and momentum are exactly conserved at all 
vertices and due to the macroscopic propagation distances of the neutrinos and muons all 
these particles follow essentially classical trajectories {i.e. corresponding to the minima 
of the classical action) which are rectilinear paths with constant velocities. The essential 
formula of Feynman's version of quantum mechanics, to be employed in the calculations 

1 The interference term is proportional to cos^i2 or sin 2 
2 Units with % = c = 1 are used throughout. 
Postulate 1 and Eqn.(7) of Ref. [12]. 

4 Heisenberg remarked that the fundamental formula (1.2) must be distinguished from that where the 
summation over intermediate states is made at the level of probabilities, rather than amplitudes, and 
that the distinction between the two formulae is 'the centre of the whole quantum theory'. 
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presented below, is [12, 15]: 

Pfi 
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(1.2) 



where P/j is the probability to observe a final state /, given an initial state i, and kj, j = 
l,n are (unobserved) intermediate quantum states. In the applications to be described 
in this paper, which, for simplicity, are limited to the case of the first two generations of 
leptons, Eqn(1.2) specialises to 5 : 



fc=l,2 



e l^fcX^fc, x D \u k , x k )(is k \n + )(iT + , x k \n + , x ) (vr + , x \S n ) 



;i.3) 



for the case of neutrino oscillations and 



E ( e+ \f I k)(f I ki X D\ni,Xk)(ni \lT + )(lT + ,X k \n + ,X )(lT + ,X \S 7r ) 
k=l,2 



(1.4) 



for the case of muon oscillations. P e -,r+ is the probability to observe the charged current 
neutrino interaction: (ui, V2)n — > e~~p following the decay: tt + — > /x + (^i, 1*2), while -P e + 7r + is 
the probability to observe the decay fi + — > e + (z/ l5 1^2) (^1, ^2)) after the same decay process. 
In Eqns.(1.3),(1.4) \u k ), k = 1, 2 are neutrino mass eigenstates while l/xjj" ), fc = 1, 2 are the 
corresponding recoil muon states from pion decay. (u k \n + ), (^ k \7t + ) and (e + \fi k ) denote 
invariant decay amplitudes, (e~\v k ) is the invariant amplitude of the charged current 
neutrino interaction, (p, x%\p, x±) is the invariant space-time propagator of particle p = 
v, n,7i between the space-time points x\ and X2 and (tt + , Xo\S n ) is an invariant amplitude 
describing the production of the n + by the source S n and its space-time propagation to the 
space-time point xo- An important feature of the amplitudes appearing in Eqns(1.3) and 
(1.4) is that they are completely defined in terms of the physical neutrino mass eigenstate 
wavefunctions \u k ). This point will be further discussed in Section 5 below. 

The difference of the approach used in the present paper to previous calculations 
presented in the literature can be seen immediately on inspection of Eqns(1.3) and (1.4). 
The initial state 6 is a pion at space time-point xo, the final state an e~ or e + produced 
at space-time point xd- These are unique points, for any given event and do not depend 
in any way on the masses of the unobserved neutrino eigenstates propagating from x k 
to xd in Eqn(1.3). On the other hand the (unobserved) space-time points x k at which 
the neutrinos and muons are produced do depend on k. Indeed, because of the different 
velocities of the propagating neutrino eigenstates, only in this case can both neutrinos 
and muons ( representing alternative classical histories of the decaying pion) both arrive 
simultanously at the unique point xd where the neutrino interaction occurs (Eqn(1.3)) or 
the muon decays (Eqn(1.4)). 

The crucial point in the above discussion is that the decaying pion, via the different 
path amplitudes in Eqns(1.3) and (1.4), interferes with itself. To modify very slightly 

5 In Eqns(1.3) and (1.4) an additional summation over unobserved states, with different physical masses 
of the decay muon, is omitted for simplicity. See Eqns.(2.1) and (2.35) below. 

6 The pion production and propagation amplitude {tt + ,xq\S- k ) contributes only a multiplicative con- 
stant to the transition probabilites. The initial state can then just as well be defined as 'pion at x \ 
rather than \S V ). This is done in the calculations presented in Section 2 below. 
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Dirac's famous statement 7 : 'Each pion then interferes only with itself. Interference be- 
tween two different pions never occurs'. 



Because of the different possible decay times of the pion in the two interfering path 
amplitudes, the pion propagators (ir + , Xk\ir + , x ) in Eqns(1.3) and (1.4) above give im- 
portant contributions to the interference phase. To the author's best knowledge, this 
effect has not been taken into account in any previously published calculation of neutrino 
oscillations. 



The results found for the oscillation phase are, for pion decays at rest: 



012 = 012 = 7 2 2W X ' 5 

(mi - m 2 ) 2 



and for pion decays in flight: 



5 12 



where 



(ml - m 2 l )E v cos9 v 

2m 2 l Am 2 (m 2 l E 7T — m^E^L 
( m 2 _ m 2) 2jB 2 cos ^ 



A 2 2 2 

Am = m 1 — m 2 



;i. 6 ) 
;i.7) 



The superscripts indicate the particles whose propagators contribute to the interference 
phase. Also E n , E v and E^ are the energies of the parent ir and the decay v and /i 
and 9 U , 9^ the angles between the pion and the neutrino, muon flight directions. In 
Eqns.(1.5) to (1.7) terms of order mf, m\, and higher, are neglected, and in Eqns.(1.6) 
and (1.7) ultrarelativistic kinematics with E n ^ ^> is assumed. Formulae for the 

oscillation phase of neutrino oscillations following muon decays or nuclear /5-decays at 
rest, calculated in a similar manner to Eqn(1.5), are given in Section 3 below. 

A brief comment is now made on the generality and the covariant nature of the cal- 
culations presented in this paper. Although the fundamental formula (1.2) is valid in 
both relativistic and non-relativistic quantum mechanics, it was developed in detail by 
Feynman [12, 14] only for the non-relativistic case. For the conditions of the calculations 
performed in the present paper (propagation of particles in free space) the invariant space- 
time propagator can either be derived (for fermions) from the Dirac equation, as originally 
done by Feynman [9] or, more generally, from the covariant Feynman path integral for 
an arbitary massive particle, as recently done in Ref . [7] . In the latter case, the invariant 
propagator for any stable particle with pole mass m, between space-time points Xi and 
Xf in free space is given by the path integral [7]: 



where r is the proper time of the particle. By splitting the integral over x(t) on the right 
side of Eqn(1.8) into the product of a series of infinitesimal amplitudes corresponding 

7 'Each photon then interferes only with itself. Interference between two different photons never oc- 
curs' [16]. 
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to small segments, At, Gaussian integration may be performed over the intermediate 
space-time points. Finally, integrating over the proper time r, the analytical form of the 
propagator is found to be [7] : 



and Hi is a first order Hankel function of the second kind, in agreement with Ref. [9]. 

In the asymptotic region where s ^> m _1 , or for the propagation of on-shell parti- 
cles [7], the Hankel function reduces to an exponential and yields the configuration space 
propagator ~ exp(— ims) of Eqn(2.11) below. It is also shown in Ref. [7] that energy and 
momentum is exactly conserved in the interactions and decays of all such 'asymptotically 
propagating' particles. The use of quasi-classical particle trajectories and the require- 
ment of exact energy-momentum conservation are crucial ingredients of the calculations 
presented below. 

The structure of the paper is as follows. In the following Section the case of neutrino 
or muon oscillations following pion decay at rest is treated. Full account is taken of the 
momentum wave-packets of the progagating neutrinos and muons resulting from the Breit- 
Wigner amplitudes describing the distributions of the physical masses of the decaying 
pion and daughter muon. The corresponding oscillation damping corrections and phase 
shifts are found to be very small, indicating that the quasi-classical (constant velocity) 
approximation used to describe the neutrino and muon trajectories is a very good one. 
The incoherent effects, of random thermal motion of the source pion, and of finite source 
and detector sizes, on the oscillation probabilities and the oscillation phases, are also 
calculated. These corrections are found to be small in typical experiments, but much 
larger than those generated by the coherent momentum wave packets. In Section 3, 
formulae are derived to describe neutrino oscillations following muon decay at rest or 
the /9-decay of radioactive nuclei. These are written down by direct analogy with those 
derived in the previous Section for pion decay at rest. In Section 4, the case of neutrino 
and muon oscillations following pion decay in flight is treated. In this case the two- 
dimensional spatial geometry of the particle trajectories must be related to the decay 
kinematics of the production process. Due to the non-applicability of the ultrarelativistic 
approximation to the kinematics of the muon in the pion rest frame, the calculation, 
although straightforward, is rather tedious and lengthy for the case of muon oscillations, 
so the details are relegated to an appendix. Finally, in Section 5, the positive aspects 
and shortcomings of previous treatments in the literature of the quantum mechanics of 
neutrino and muon oscillations are discussed in comparison with the method and results 
of the present paper. 

8 The metric for four-vector products is time-like. 




(1.9) 



where 8 : 
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2 Neutrino and muon oscillations following pion de- 
cay at rest 



To understand clearly the different physical hypotheses and approximations underly- 
ing the calculation of the particle oscillation effects it is convenient to analyse a precise 
experiment. This ideal experiment is, however, very similar to LNSD [17] and KAR- 
MEN [18] except that neutrinos are produced from pion, rather than muon, decay at 
rest. 

The different space-time events that must be considered in order to construct the 
probability amplitudes for the case of neutrino oscillations following pion decay at rest 
are shown in Fig.l. A n + passes through the counter Ca, where the time t is recorded, 
and comes to rest in a thin stopping target T (Fig. la)). For simplicity, the case of only two 
neutrino mass eigenstates v\ and V2 of masses m\ and rri2 {mi > m^) is considered. The 
pion at rest constitutes the initial state of the quantum mechanical probability amplitudes. 
The final state is an e~p system produced, at time to, via the process (vi, v-i)n — > e~p at a 
distance L from the decaying rt + (Fig. Id)). Two different physical processes may produce 
the observed e~p final state, as shown in Fig. lb) and lc), where the pion decays either 
at time t\ into v\ or at time £2 into V2- The probability amplitudes for these processes 
are, up to an arbitary normalisation constant, and neglecting solid angle factors in the 
propagators: 



Note that following the conventional '/i' (final, initial) ordering of the indices of ma- 
trix elements in quantum mechanics, the path amplitude is written from right to left 
in order of increasing time. This ensures also correct matching of 'bra' and 'kef sym- 
bols in the amplitudes. In Eqn(2.1), < e~p\T R \nv i >, < z/j / u + |Tft|7r + > are 'reduced' 
invariant amplitudes of the v charged current scattering and pion decay processes, re- 
spectively, BW(W fl ( i ),m^ Tfj) and BW(W 7V: m 7V , T w ) are relativistic Breit-Wigner ampli- 
tudes, U ei and U ifl are elements of the unitary Maki-Nagagawa-Sakata (MNS) matrix [19], 
U a i, describing the charged current coupling of a charged lepton, a, ( a — e,fi,r) to 
the neutrino mass eigenstate i. The reduced invariant amplitudes are defined by fac- 
toring out the MNS matrix element from the amplitude for the process. For example: 
< e~p\T\ni>i >= U ei < e~p\T R \ni>i >. Since the purely kinematical effects of the non- 
vanishing neutrino masses are expected to be very small, the reduced matrix elements may 
be assumed to be lepton flavour independent: < e~p\T R \nvi >~< e~p\T R \nuo > where uq 
denotes a massless neutrino. In Eqn(2.1), D is the Lorentz-invariant configuration space 
propagator [7, 9] of the pion or neutrino. The pole masses and total decay widths of the 
pion and muon are denoted by m n , T n and m M , T M respectively. For simplicity, phase 
space factors accounting for different observed final states are omitted in Eqn(2.1) and 
subsequent formulae. 

Because the amplitudes and propagators in Eqn(2.1) are calculated using relativistic 
quantum field theory, and the neutrinos propagate over macroscopic distances, it is a 




< e p\T R \nVi > U ei D(x f - x { , t D - t h m^BWiW^.m^ Y ^) 



xBWiW^m^Y^dW^ 2 = 1,2 



(2.1) 
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good approximation, as already discussed in the previous Section, to assume exact energy- 
momentum conservation in the pion decay process, and that the neutrinos are on their 
mass shells, i.e. p 2 = mf, where Pi is the neutrino energy-momentum four-vector. In 
these circumstances the neutrino propagators correspond to classical, rectilinear, particle 
trajectories. 

The pion and muon are unstable particles whose physical masses W n and W^m differ 
from the pole masses m n and m M appearing in the Breit-Wigner amplitudes and covariant 
space-time propagators in Eqn(2.1). The neutrino momentum Pj will depend on these 
physical masses according to the relation: 



P 



[W 2 _ ( m . + w^) 2 ][w 2 - (m - w^f] 

2W n 



(2.2) 



Note that, because the initial state pion is the same in the two path amplitudes in 
Eqn(2.1) W„ does not depend on the neutrino mass index %. However, since the pion decays 
resulting in the production of v\ and v-i are independent physical processes, the physical 
masses of the unobserved muons, W^, i = 1,2, recoiling against the two neutrino 
mass eigenstates are not, in general, the same. In the following kinematical calculations 
sufficient accuracy is achieved by retaining only quadratic terms in the neutrino masses, 
and terms linear in the small quantities : S n = — m n , 5i = — m^. This allows 

simplification of the relativistic Breit-Wigner amplitudes: 



BW(W,m,T) = 



Tm 



W 2 — m 2 + imT 
Tm 

5(2m + 5)+ imT 



2(5 + 
BW(5,T) + 0(5 2 



(2.3) 



Developing Eqn(2.2) up to first order in m 2 , 5i and 5 n yields the relation : 



P 



1 - 



m 2 (ml + m 2 ) 5 W (m 2 n + m 2 ) 



+ 



(ml — m 2 ) 2 m n (ml 



m. 



25jm M 5 n m 2 (ml + m 



+ 



ml -m 2 ^ m n (ml-m 2 i y 



(2.4) 



where 



Pn = 



ml -ml 



2m Tl 



29.8MeV 



(2.5) 



The term ~ S^m 2 which is also included in Eqn(2.4) gives a negligible 0(mf) contribution 
to the neutrino oscillation formula. For muon oscillations, however, it gives a term of 
O(mf) in the interference term, as discussed below. Similarly, the exact formula for the 
neutrino energy: 



Wi - Wfa + m. 



2W„ 



(2.6) 
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in combination with Eqn(2.4) gives, for the neutrino velocity: 



Pi m 



2 



Ei 2P 2 



7-771,, 



+ OtwZXJi) (2-7) 



m^(m\ - m 2 ^ m^-m^ 

This formula will be used below to calculate the neutrino times-of- flight: t{ 1 % — 1,2. 

From the unitarity of the MNS matrix, the elements C/ a j may be expressed in terms 
of a single real angular parameter 9: 

U el = U le = U ia = U 2li = cose (2.8) 
U e2 = U 2e = -U lil = -U lli = wn9 (2.9) 

The parts of the amplitudes requiring the most careful discussion are the invariant 
space-time propagators D, as it is mainly their treatment that leads to the different 
result for the neutrino oscillation phase found in the present paper, as compared to those 
having previously appeared in the literature. In the limit of large time-like separations, 
the propagator may be written as [7, 9]: 



m 



D(Ax,At,m) = expHm v /(At)2 - (Ax)*} (2.10) 

2mJ(At) 2 -(Ax) 2 i 



D is the amplitude for a particle, originally at a space-time point (x*i, ti), to be found at 
(xf, tf) and Ax = x/ — Xj, At = tf — U. In the following, according to the geometry of the 
experiment shown in Fig.l, only one spatial coordinate will be considered (Ax = Xf — Xj) 
and only the exponential factor in Eqn(2.10), containing the essential phase information 
for particle oscillations will be retained in the amplitudes. Solid angle correction factors, 
taken into account by the factor in large brackets in Eqn(2.10), are here neglected, but 
are easily included in the final oscillation formulae. Writing then 



D(Ax, At, m) ~ exp[-im^J (At) 2 - (Ax) 2 } = exp[-imAr] = exp[-iA0] (2.11) 

it can be seen that the increment in phase of the propagator, A0, when the particle 
undergoes the space-time displacement (Ax, At) is a Lorentz invariant quantity equal 
to the product of the particle mass and the increment, At, of proper time. Using the 
relativistic time dilatation formula: 

E 

At = 7 Ar = — At (2.12) 
m 

and also the relation, corresponding to a classical, rectilinear, particle trajectory: 

At = — = —L (2.13) 

v p 
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gives, for the phase increments corresponding to the paths of the neutrinos and the 
pion in Fig.l: 



2 2 
IT) TO 

Atf = m^n = -f At, = -±L 



1 - 



5 W (ml + mj) 2^m 



+ 



(m^ - 777^) ml - m 2 



(2.14) 



m^U - t ) = m n (t D - t ) 
m w (t D - t ) - m^L \ 1 



777J 



2P 2 



45j777 M 



28* {ml + ml) 
m n (ml — m 2 ^) ' ml — m 2 ^ 



(2.15) 



where terms of 0(t77^) and higher are neglected. 



Making the substitution t; L — 1 — > t D — t — L/vi in the exponential damping factor due 
to the pion lifetime in Eqn(2.1) and using Eqns(2.3),(2.11),(2.14) and (2.15) Eqns(2.1) 
may be written as: 



A,; 



I 



< e p\T R \nu > U ei U ilM < u fi + \T R \n + > 



2 (ft + 



YJs^+W)' 



3 i<fo-^(t D -to-tf) 



exp i 



m~ / m 



Po V2P 



1 L 



d5i 



where 



on = 



0o = m w (L -t D + t ) 
Amjm^m 7T (ml + m 2 )L 
(ml — m 2 ^) 3 
2m 2 (ml + mD 2 L 



a n (i) 



(ml - mlf 



and 



tf^(i + J) + oW) 



m; 



1,2 (2.16) 

(2.17) 
(2.18) 

(2.19) 
(2.20) 



In Eqns(2.18) and (2.19) imaginary parts of relative size ~ T^/m* ~ 2.0 x 10 16 are 
neglected. 



To perform the integral over Si in Eqn(2.16) it is convenient to approximate the mod- 
ulus squared of the Breit-Wigner amplitude by a Gaussian, via the substitution: 



r 

2 



S + i\ 



r / jj^ 
2 i^T? 



2 r 

-((J-i-)exp 



(2.21) 



where the width of the Gaussian is chosen so that it has approximately the same full 
width at half maximum, T, as the Breit-Wigner function. After the substitution (2.21), 
the integral over Si in Eqn(2.16) is easily evaluated by a change of variable to 'complete 
the square' in the argument of the exponential, with the result: 



2 

f: 



: ))exp 



3§ 



+ ia.iSi dSi 



r M exp 



12 



otiT 



- 1 
(2.22) 
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Eqn(2.18) gives, for ctj, the numerical value: 

a, = 3.1 x 10 14 ( — V L(m) MeV" 1 (2.23) 

For typical physically interesting values (see below) of — 1 eV and L = 30 m, on takes 
the value 0.48 MeV -1 , so that 

aiT^ = 0.48 x 3.00 x 10~ 16 = 1.4 x 10~ 16 



Then, to very good accuracy, I± — I 2 — —iyiv/3T^, independently of the neutrino mass. 
It follows that for neutrino oscillations, the muon mass dependence of the amplitudes may 
be neglected for any physically interesting values of rrii and L. 

From Eqns (2.16) and (2.22) the probability to observe the reactions (vi, v-i)n — > e~p 
at distance L from the pion decay point and at time to is: 

P(e-p\L,t D ) = \A 1 + A 2 \ 2 
ttF 2 

= -j 1 ] < e~p\T R \nv > | 2 | < v /i + \T R \ir + > \ 2 

r 2 

xsin 2 flcos 2 fle- r ^ D ~' o) ^ (2.24) 

4 + 



x <^ e * 1 + e n 2 



-2e n ^ i? e exp i 



Am 2 



2P n 



1 L a n (2)}S n 



The time dependent exponential factors in the curly brackets of Eqn(2.24) are easily 
understood. If m\ > ra 2 then t{ 1 > t 2 l . This implies that the neutrino of mass mi 
results from an earlier decay than the neutrino of mass m 2 , in order to be detected at the 
same time. Because of the exponential decrease with time of the pion decay amplitude, 
the contribution to the probability of the squared amplitude for the neutrino of mass 
mi is larger. The interference term resulting from the product of the decay amplitudes 
of the two neutrinos of different mass, has an exponential factor that is the harmonic 
mean of those of the squared amplitudes for each neutrino mass eigenstate, and so is also 
suppressed relative to the squared amplitude for the neutrino of mass m\. The integral 
over the physical pion mass is readily performed by replacing the Breit-Wigner function 
by a Gaussian as in Eqn(2.21). This leads to an overall multiplicative constant ^-K/2>T n 
and a factor: 

F V (W V ) = exp[-K(l) - a w (2)) 2 r 2 /12 (2.25) 

multiplying the interference term. For Am 2 = (leV) 2 and L = 30m the numerical value 
of this factor is exp(— 1.3 x 10~ 29 ). This tiny correction is neglected in the following 
equations. 

Integrating over to gives the average probability to observe the e~p final state at 
distance L: 

P(e~p\L) = ^T7#l < e-p\T R \nu > | 2 | < iy fi + \T R \n + > | 2 sin 2 #cos 2 # 

r,m 2 Am 2 2m n m 2 Am 2 \ 
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Where all kinematical quantities are expressed in terms of Am 2 , m n and m M . Note that 
the minimum value of to is t + t{ 1 , t + t{ 1 and t +t[ l for the squared amplitude terms of 
neutrinos of mass m±, m 2 and the interference term, respectively. On integrating over to, 
the squared amplitude terms give equal contributions, the larger amplitude for mass mi 
being exactly compensated by a smaller range of integration. The exponential damping 
factor in the interference term in Eqn(2.26) is derived using the relations: 

t{ 1 -tfl=(-L-- J—) L ~ (v 2 (u) - Vl {v))L (2.27) 



and 



to obtain 



2 

mi 



M = l-515 + OK) < = 1>2 (2.28) 



2P 2 



« « _ (mj - m 2 )L _ Am 2 L 

The damping factor arises because the difference in the times-of-flight of the two neutrino 
paths is limited by the pion lifetime. It will be seen below, however, that for distances L 
of practical interest for the observation of neutrino oscillations, the damping effect is tiny. 

The part of the oscillation phase in Eqn(2.24) originating from the neutrino propa- 
gators (the term associated with the '1' within the large curved brackets) differs by a 
factor two from the corresponding expression in the standard formula. The contribution 
to the oscillation phase of the propagator of the decaying pion ( the term associated with 
m 7r /(2P ) in the large curved brackets of Eqn(2.24)) has not been taken into account in 
any published calculation known to the author of the present paper. The oscillation phase 
in Eqn(2.26) is 2m 2 / (m 2 — m 2 ) = 2.685 times larger than that given by the standard 
formula (1.1). For L = 30m, as in the LNSD experiment, the first oscillation maximum 
occurs for Am 2 = 0A6(eV) 2 . Denoting by 4>i2 the phase of the cosine interference term 
in Eqn(2.26), the pion lifetime damping factor can be written as: 

F»(r n ) = exp (-F^^ = exp(-1.58 x lO" 18 ^) (2-30) 

so the damping effect is vanishingly small when (fy^ — 1. 

The oscillation formula (2.26) is calculated on the assumption that the decaying pion 
is at rest at the precisely defined position x«. In fact, the positive pion does not bind with 
the atoms of the target, but will rather undergo random thermal motion. This has three 
effects: an uncertainy in the value of Xi, a Doppler shift of the neutrino energy and a time 
dilatation correction correction factor of 1/7^ in the equation (2.15) for the pion phase 
increment. Assuming that the target is at room temperature (T= 270° K), the mean 
kinetic energy of 3kT/2 correponds to a mean pion momentum of 2.6 x 10~ 3 MeV and 
a mean velocity of ~ 5.6 km/sec. The pion will move, in one mean lifetime (2.6 x 10~ 8 
sec), a distance of 146/xm. This is negligible as compared to L (typically > 30m) and so 
Eqn(2.26) requires no modification to account for this effect. 

The correction factor due to the Doppler effect and time dilatation is readily calculated 
on the assumption of a Maxwell-Boltzmann distribution of the pion momentum: 
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Here p n = ^/2kTm n = 2.64 x 10" 3 MeV. Details of the calculation are given in Appendix 
A. The interference term in Eqn(2.26) is modified by a damping factor: 



F u (Dop) = exp 




p 7r Am 2 \m 



2m ff P LP C 



o 



'7T 



1 L 



) 



2 



I 



(2.32) 



while the argument of the cosine term acquires an additional phase factor: 



F(Dop) 



3 / p w \ 2 Am 2 r 3m. 



- 1 L 



(2.33) 



4 \mj P 12P 



For 4$ = 1, F u (Dop) = 1 - 6.7 x 1CT 10 and ^(Dop) = 1.2 x 1CT 9 . 

If the target in which the pion stops is of thickness £t, then the effect of different 
stopping points of the 7r (assumed uniformly distributed) is to multiply the interference 
term in (2.26) by the factor: 



If the position of the neutrino interaction point within the target has an uncertainy of 
±£d/2 a similar correction factor is found, with the replacement It — > £d in Eqn(2.34). 
The calculation of this correction factor is also described in Appendix A. 

If the target T in which the pion comes to rest (Fig. la)) is chosen to be sufficiently thin, 
the pion decay process may be detected by observing the recoil muon in the counter Ca 
at times t± or t 2 (Fig. lb) or lc)). A sufficiently accurate measurement of the times ti, t 2 
and to would, in principle, enable separation of the different processes in Figs, lb) and lc) 
by the observation of separated peaks in the time-of-flight distribution at t{ = t D — 1\ and 
t{ = t D — t 2 . In this case the interference term in Eqn(2.26) vanishes as the two alternative 
space-time paths leading to the neutrino interaction are distinguishable. However, for 
L = 30m and Am 2 = l(eV) 2 the difference in the times of flight is only 5.6 x 1CT 23 
sec, more than ten orders of magnitude smaller than can be measured with existing 
techniques. As discussed in Section 5 below, the momentum smearing due to the Doppler 
effect at room temperature is some eleven orders of magnitude larger than the shift due 
to a neutrino mass difference with Am 2 = l(eV) 2 . Thus, even with infinitely good time 
resolution, separation of such neutrino mass eigenstates by time-of-flight is not possible. 

The ideal experiment, described above to study neutrino oscillations, is easily adapted 
to the case of oscillations in the decay probability of muons produced by charged pion 
decay at rest. As previously pointed out in Ref. [8], such oscillations will occur if neutrinos 
with different masses exist. As before, the pion stops in the target T at time to ( see 
Fig. 2a)). At time t\ the pion decays into v\ and the corresponding recoil muon (/xi), whose 
passage is recorded in the counter Cb (Fig. 2b)). Similarly, a decay into u 2 and fi 2 may 
occur at time t 2 (Fig. 2c)). With a suitable choice of the times ti and t 2) such that muons 
following the alternative paths both arrive at the same time t D dX the point Xf, interference 
occurs between the path amplitudes when muon decay occurs at the space-time point (xf, 
to) in the detector D (Fig. 2d)). The probability for two classical trajectories to arrive 
at exactly the same space-time point of course vanishes. The correct way to consider the 
quantum mechanical calculation is rather to ask given that the muon decay occurs at the 




(2.34) 
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point {x, to), does the muon recoil against v\ or z/ 2 ? If these two possiblities are not 
distinguished by the measurement of the decay process, the corresponding probability 
amplitudes (not probabilities) must be added in the calculation of the probability of the 
observed decay process. The path amplitudes corresponding to muons recoiling against 
neutrinos of mass m\ and m 2 are: 



(kl)i 



< e + u k ui\T\fx + > e 27 » D(xf — Xi,t D 
BW(W^,m^r^U^ < v^ + \T r \-k + > 
e -¥(* i -*o) jD (o, ti - t , m 7V )BW(W 7V , m w , T^dW^ i = 1, 2 (2.35) 

The various factors in these equations are defined, mutatis mutandis, as in Eqn(2.1). 



With the same approximations, concerning the neutrino masses and the physical pion 
and muon masses, as those made above, the velocity of the muon recoiling against the 
neutrino mass eigenstate Vi is: 



1 - 



[ml — m 2 ^) 2 {ml + mV) 



A5 i m II ml 
mi -ml 



d^mfmlm^Sm 2 ^ - ml) 
W - ^) 2 (ml - ml) 



(2.36) 



where 



ml 



m. 



ml + m 2 ^ 



(2.37) 



Comparing with Eqn(2.7), it can be seen that for the muon case, unlike that where 
neutrino interactions are observed, there are pion and muon mass dependent correction 
terms that are independent of the neutrino masses, implying a velocity smearing effect 
due to the physical pion and muon masses that is ~ ml/m 2 larger than for the case of 
neutrino oscillations. 

The phase increments corresponding to the paths of the muons and the pion in Fig. 2 
are, using (2.4) 9 and (2.12)-(2.15) and (2.36): 



A0f = 



m 2 L 

fjj 



m 2 L 



1 + 



Pt Po 
5 n m 2 (ml + ml) 



2m n P 2 



5 W 



{ml + ml) 



m w (ml — m 2 u 



m Tl 



{ml - ml) 2 



2(^m M 
ml — m 2 



(2.38) 



A0I 



m*(t D - 1 ) - 



m w {ti - t ) 
m n L 



m n {t D - 1 ) - 



m w L 



1 + 



Amfmlml 



{ml -ml) 2 {ml + ml) 



AS^m^ml ASjin^ml 8<J w 77i?7nJmJ(3mJ - ml) 



m%-m* 



m^ -ml 



{mi - ml) 2 {ml - m*) 



3 Note that, in the pion rest frame Pi = Pj*. 



(2.39) 
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where 



ES 



ml + ml 
2m n 



(2.40) 



Using Eqns(2.11),(2.38) and (2.39) to re- write the space-time propagators in Eqn(2.35), 
as well as Eqn(2.3) for the Breit-Wigner amplitudes gives: 



A 



(ki)i 



J < e + u k u l \T\/2 + > e ~^~U^ < u^ + \T\ti + > 



2 (Si + i^) 2 (8* + itf 



2P| 



1 - 



where 



, u 



-tn + t ) - 



Am^m^L 
ml -ml 



1,2 



2m 2 L 



m2 - m 2 



m\{hm% - \\m\m\ - m^m^ - mj) 



K-^)K-^) : 



and 



>(i) 



= 14 



mfmy 
2m 7T P§ i 



(2.41) 

(2.42) 
(2.43) 

(2.44) 
(2.45) 



where, as in Eqns(2.18) and (2.19), imaginary parts of order T^/m^ are neglected. Making 
the substitution (2.21) and performing the integral over Si according to Eqn(2.22), the 
following formula is found for the probability for muon decay at distance L and time to- 



P(e + u k ui\L,t D ) 



\A 



(/') 



;.y)l + \kl)2 



[1 - < e + v k V t \T\n + > 



\ 2 e < 



< u /j, + \T R \7r + > | 2 e - r -(^-*<») 



m + T 4) 



jsin 2 Oe^V) + cos 2 9e 
-2 sin 9 cos e ^ ( *£(i) +t £(a) 
i?e exp i 



r ^ 



(2.46) 



Where the effect of the non-zero neutrino masses are neglected in the reduced pion decay 
amplitude so that < i/j// + |7jj|7r + >~< uofi + \Tf>\ii + > and this amplitude is a common 
factor in both path amplitudes. The muon path difference yields the term associated with 
Eq /m^ in the interference phase in Eqn(2.46) while the pion path is associated with '1' 
in the large round brackets. The numerical value of the damping factor: 



F^) = exp[-^^][l-^i] 2 



(2.47) 
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resulting from the integral over the physical muon mass is, for L = 30m, 0.774, so, unlike 
for the case of neutrino oscillations, the correction is by no means negligible. This is 
because, in the muon oscillation case, the leading term of is not proportional to the 
neutrino mass squared. The non-leading terms proportional to m 2 have been neglected 
in Eqn(2.43). This correction however effects only the overall normalisation of the oscil- 
lation formula, not the functional dependence on L arising from the interference term. 
Integrating over 5 W using Eqns(2.21) and (2.22), as well as over to, gives the probability of 
muon decay, into the final state e + i> k Vi, at distance, L, from the production point, where 
all kinematical quantities are expressed in terms of Am 2 , m n and m M : 



P(e + ^|L) 




4 T^m^m w 2 

1 ~~ O 2 2^ 

3 m A — m z 



2T ll m 7T m fl (ml - m 2 



{ml + m 2 ) 3 



2T n m 2 m 2 Am 2 2m n m 2 Am 2 

sin2#exp[ — L] cos - — ^ — - — tttitL 



[mi 



m 



\mt 



m 



2\2 




[A/ \ 7T fl/ 

In this expression the correction due to the damping factor of the interference term: 



= exp[-K(l) - <(2)) 2 r 2 /12] 



(2.49) 



arising from the integral over the physical pion mass has been neglected. For Am 2 = 
(leV) 2 and L = 30m the numerical value of this factor is exp(— 5.2 x 10 -30 ). Denoting by 
012^ the argument of the cosine in Eqn(2.48), the exponential damping factor due to the 
pion lifetime may be written as: 



F»(r n ) = exp 




(2.50) 



For <pi2 — 1) F^iY^) = exp[— 4.4 x 10 ] so, as in the neutrino oscillation case, the pion 
lifetime damping of the interference term is very small. 



Introducing the reduced muon decay amplitude: 
< e + u k T7 l \T\fi + >= U ek U^i < e + u k T7 l \T R \fj J + >~ U ek U^ < e + v V \T R \n + > 



(2.51) 



the total muon decay probability is given by the incoherent sum over the four possible 
final states containing massive neutrinos: 

P(e + uV\L) = ]T^P(e + ^|L) 

k=l 1=1 

= jz\Uek\ 2 jz\U,i\ 2 P{e + v Q V \L) 

k=l 1=1 

= P(e + u V \L) (2.52) 

where the unitarity of the MNS matrix has been used. P(e + z/ !7o|Z-) is given by the 
replacement of < e + u k U l \T\fi + > by < e + u U \T R \fi + > in Eqn(2.48). Eqn(2.52) shows 
that the muon decay width is independent of the values of the MNS matrix elements. 
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Corrections due to time dilatation and the Doppler effect are calculated in a similar 
way to the neutrino oscillation case with the results (see Appendix A): 



F^{Dop) = exp 



' p 7T mf l Am 2 VQ 



"3 








')') 


2 


m 7T _ 





and 



3 


(2jL 


* 2 m^Am 2 


f ^ 


2 


Vm,, 




2m,,- 



(2.53) 



(2.54) 



As for neutrino oscillations, the corresponding corrections are very small for oscillation 
phases of order unity. 

The phase of the cosine in the interference term is the same in neutrino and muon 
oscillations, as can be seen by comparing Eqns(2.26) and (2.48). It follows that the target 
or detector size correction (Eqn(2.34)) is the same in both cases. 

Neutrino and muon oscillations from pion decay at rest then have an identical oscil- 
lation phase for given values of Am 2 and L. In view of the much larger event rate that 
is possible, it is clearly very advantageous in this case to observe muons rather than neu- 
trinos, since the rate of neutrino oscillation events is severely limited by the very small 
neutrino interaction cross section. In fact, it it not necessary to observe muon decay, as 
in the example discussed above. The oscillation formula applies equally well if the muons 
are observed 10 at the distance L using any high efficiency detector such as, for example, 
a scintillation counter. According to Eqn(1.2), interference between the path amplitudes 
must occur if the muon detection device does not discriminate muons recoiling against v\ 
from those recoiling against vi- 



3 Neutrino oscillations following muon decay or beta- 
decay at rest 

The formula describing — > V e neutrino oscillations' 11 following the decay at rest of a 
ix + — > e + (z/x, U2){Vi,V2) is easily derived from the similar formula for n + decay at rest, 
(2.25). Because the neutrino momentum spectrum is continous, smearing effects due to 
the finite muon lifetime may be neglected from the outset. The phase increment associated 

10 Note that, in this case, the final state of the path amplitude is that of the detection process by 
which the muon is recorded. As is the neutrino in neutrino oscillations, the muon itself contributes an 
unobserved intermediate state in the general path amplitude formula (1.2) 

n This experiment is also conventionally termed '77 e appearence'. As discussed in more detail in Section 
5 below, '77 e ' and do not exist, as physical states, if neutrinos are massive and the MNS matrix is 
non-diagonal. It is still, however, current practice in the literature to use the symbols V e ' , l i>^ and 
'z/ r ' to refer to massive neutrinos. This is still a useful and meaningful procedure if it is employed only 
to identify, in a concise manner, the type of charged current interaction by which the neutrinos are 
produced or detected, i.e. 'zV means neutrinos (actually several, with different generation numbers) 
produced together with the charged lepton I or detected by observation of the charged lepton I. It should 
not be forgotten however that only the wavefunctions of the mass eigenstates vi occur in the amplitudes 
of Standard Model processes. 
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with the neutrino path is then given by Eqn(2.14) with the replacements Po — > pj and 
0, where pj is the antineutrino momentum. The phase increment of the decaying 
muon is given by the same replacements in Eqn(2.15) with, in addition, m n — > and 
— > T^. The formula, analogous to Eqn(2.26), for the time-averaged probability to 
detect the process (T7i,T7 2 )p — > e + n at a distance L from the muon decay point is then: 



P(e + n,/i\L) 



< e + n\T R \pu > | 2 | < v T7 e + \T R \n + > 



-2 sin 2 9 cos 2 9 



L r r M Am 2 

x a - exp[ ^2-^] cos 



Am 2 / m 



Pf V2P f 



1 L 



(3.1) 



The standard neutrino oscillation formula, hitherto used in the analysis of all exper- 
iments, has instead the expression Am 2 L/{2P v ) for the argument of the cosine term in 
Eqn(3.1). Denoting my Am| the value of Am 2 obtained using the standard formula, and 
Am 2 FP that obtained using the Feynman Path (FP) formula (3.1) then: 

/\ 771 ^ 

p- 

For a typical value of P v of 45 MeV, Eqn(3.2) implies that Am 2 FP ~ 2.9Am|. Thus 
the oscillation signal from fi + decays at rest reported by the LNSD Collaboration [17] 
corresponding to Am| ~ 0.5 (eV) 2 for sin 2 29 ~ 0.02 implies Am 2 FP ~ 1.5 (eV) 2 for a 
similar mixing angle. 

For the case of /9-decay: 

N(A,Z) ^ N(A,Z + I)e-(V U V 2 ) 

m n in the first line of Eqn(2.15) is replaced by Ep, the total energy release in the /9-decay 
process: 

E p = M N (A, Z) - M N (A, Z + 1) (3.3) 

where M^(A, Z) and M^(A, Z+l) are the masses of the parent and daughter nuclei. That 
the phase advance of an unstable state, over a time, At, is given by exp(— iE*At) where E* 
is the excitation energy of the state, is readily shown by the application of time-dependent 
perturbation theory to the Schrodinger equation [20] 12 . A more intuitive derivation of this 
result has also been given in Ref. [21]. In the present case, E* = Ep. Omitting the lifetime 
damping correction, which is about eight orders of magnitude smaller than for pion decay, 
given a typical /9-decay lifetime of a few seconds, the time-averaged probabilty to detect 
Vi,J7 2 via the process (Vi,V 2 )p — > e + n, at distance L from the decay point is given by the 
formula, derived in a similar way to Eqns(2.26) and (3.1): 

p , + „ in | < e + n\T R \pu > | 2 | < e -V N(A,Z + l)\T R \N(A,Z) > | 2 
P{e n, p\L) = 



x < sin 4 9 + cos 4 9 + 2 sin 2 9 cos 2 9 cos 



Am 2 / E 



(*-l 



Pu V2P F 



(3.4) 



2 See also Eqn(20) of Chapter V of Ref[16] 
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where T^ 1 = the lifetime of the unstable nucleus N(A, Z). Until now, all experiments 
have used the standard expression Am 2 L/(2P 17 ) for the neutrino oscillation phase. The 
values of Am 2 found should be scaled by the factor {Ep/P v — 2) _1 , suitably averaged over 
Ptj, to obtain the Am 2 given by the Feynman Path formula (3.4). 



Neutrino and muon oscillations following pion de- 
cay in flight 



In this Section, the decays in flight of a 7r + beam with energy E n ^> m 7r into /U + (z/i, i/ 2 ) 
are considered. As the analysis of the effects of the physical pion and muon masses have 
been shown above to give negligible corrections to the L dependence of the oscillation 
formulae, for the case of decays at rest, such effects will be neglected in this discussion of 
in-flight decays. The overall structure of the path amplitudes for neutrinos and muons is 
the same as for decays at rest (see Eqns(2.1) and (2.35)). However, for in-flight decays, in 
order to calculate the interfering paths originating at different and terminating at common 
space-time points, the two-dimensional spatial geometry of the problem must be properly 
taken into account. 

In Fig. 3 a pion decays at A into the 1 mass eigenstate, the neutrino being emitted 
at an angle 61 in the lab system relative to the pion flight direction. If m\ > m 2 a later 
pion decay into the 2 mass eigenstate at the angle 6 X + 5 may give a path such that both 
eigenstates arrive at the point B at the same time. A neutrino interaction (ui, v 2 )n — > e~p 
occuring at this space-time point will than be sensitive to interference between amplitudes 
corresponding to the paths AB and ACB. The geometry of the triangle ABC and the 
condition that the 1 and 2 neutrino mass eigenstates arrive at B at the same time gives 
the following condition on their velocities: 



V\{v) sin#2 vi{y) sin(0 2 — #i) 



v 2 {v) sin^i 



sin#i 



(4.1) 



Expanding to first order in the small quantity 5 = 9 2 — 6±, rearranging, and neglecting 
terms of 0(mf), gives: 



v 2 (u) - vi(v) 



Am 2 



2E 2 sin 6 1 



1 — V n COS 61 



where the relation: 



Vi{y) = 1 



m 2 

m 



+ 0(mt) 



has been used. Rearranging Eqn(4.2): 

Am 2 



2E 2 



v n sin 0\ 
1 — v n cos 61 



(4.2) 



(4.3) 



(4.4) 



The difference in phase of the neutrino paths AB and CB is (see Eqn(2.14)): 

m\AB m 2 CB 4 4 



Pi 



(4.5) 
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Since the angle 5 is ~ Am 2 , the difference between AB and CB is of the same order, and 
so: 



012 = 



Am 2 L 



+ 0(rri,rr4) 



(4.6) 



cos 9±E U 

where P\ ~ P 2 — E u , the measured neutrino energy. From the geometry of the triangle 
ABC: 



AC 



AB 



L 



sin 5 sin 2 cos0i sin(6 l 1 + 5) 
So, to first order in 5, and using Eqn(4.4): 

L5 Am 2 L 



and 



AC = Ax n 



At n 



cos 9 1 sin 0! 2E 2 cos0i (1 — tv cos0i) 
Ax^ Am 2 L 1 



1^ 2i?2 COS 6*i (1 — V w COS 0i) 

Eqns(4.8) and (4.9) give, for the phase increment of the pion path: 

Am 2 E 7T L (1-vl) 



A(f) n = m 7r (r 2 — Ti)— m^Ar = E^At^ - p^Ax^ = 



2E 2 cos 61 (1 — v n cos 9 1 



(4.7) 

(4.8) 
(4.9) 

(4.10) 



In Eqn(4.10), the Lorentz invariant character of the propagator phase is used. Setting 
cos 0i = 1 and v n = 0, gives for Acjf a prediction consistent with that obtained from 
Eqn(2.15). Eqns(4.6) and (4.10) give, for the total phase difference of the paths AB, 
ACB: 



I ZA7T 

>12 



A<j) AB - A<f> 



ACB 



0i 2 - 



Am 2 L 
cos9\E u 



1 - 



E, (1 - vl) 
2E U (1 -v v cos0i) 



Using the expressions, valid in the ultra-relativistic (UR) limit where v n ~ 1 : 



1 — V n COS 0i = 



m; 



1 



E 2 (1 + COS0,! 



and 



^(m 2 -m 2 ) 
" = 2m 2 ( + C ° S u) 



(4.11) 



(4.12) 



(4.13) 



where 0* is the angle between the directions of the pion and neutrino momentum vectors 
in the pion rest frame, Eqn(4.11) may be rewritten as: 

-L (4.14) 



Am 2 



m. 



cos 9 X E V (m 2 - m 2 ) 



For 0! = the oscillation phase is the same as for pion decay at rest (see Eqn(2.26)) 
since in the latter case, E v ~ P — { m W ~ m jD/(2 m 7r)- Using Eqn(4.14), the probability 
to observe a neutrino interaction, at point B, produced by the decay product of a pion 
decay occuring within a region of length, I Deo (<C L) centered at the point A, in a beam 
of energy E n , is given by a formula analagous to Eqn(2.26): 



P(e-p\L,9 1 ) 



E w 

X < 1 - COS 



< e p\T R \nu > \ 2 \ < h> ii + \T R \ix + > | 2 sin 2 0cos 2 



m 2 Am 2 L 



(mj — mfyE v cos #1 



(4.15) 
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As in the case of pion decay at rest, Eqn(2.26), the oscillation phase differs by the factor 
2m^/(mJ — mfy = 2.685 from that given by the standard formula. 

The derivation of the formula describing muon oscillations following pion decays in 
flight is very similar to that just given for neutrino oscillations. The condition on the 
velocities so that the muons recoiling against the different neutrino mass eigenstates arrive 
at the point B (see Fig. 3) at the same time, is given by a formula analagous to (4.2): 



M/i)M/i) -^cosfli] 

Av(ii) = v 2 (n) - ui(fi) = — 6 

v n sin t>i 



(4.16) 



The formula relating the neutrino masses to the muon velocities is, however, more difficult 
to derive than the corresponding relation for neutrinos, (4.3), as the decay muons are 
not ultra-relativistic in the pion rest frame. The details of this calculation are given in 
Appendix B. The result is: 



Av(fj,) 



mz(m 2 + m^)Am 
E 2 (ml -ml) 2 



2m 2 t E 7r 



(ml + m 2 l )E IXi 



+ 0(mi,m 4 2 ) 



Using Eqn(4.16) and the relation, valid to first order in 5,: 

Lb 

At = 



v v cos 9\ sin Oi 



(4.17) 



(4.18) 



where At is the flight time of the pion from A to C in Fig. 3 (and also the difference in 
the times-of- flight of the muons recoiling against the two neutrino eigenstates), the angle 
8 may be eliminated to yield: 



At 



Av{p)L 



V\ (n) cos 9\ [vi (//) — v n cos 6i 
Using now the kinematical relation (see Appendix B): 

(ml 



Vi(jj) — V n COS 6>! = 



m. 



{ml + m 2 )E u 



and the expression for the phase difference of the paths AB and ACB: 

A(f) AB - A(f) ACB 



At - — N 



E 



(4.19) 



(4.20) 



(4.21) 



together with Eqn(4.19), it is found, taking the UR limit, where fi(/i),iv ~ 1, that 

(mlE n - mlE,)L 



'12 



2m 2 Am 2 

fx 



E 2 {m 2 n - m 2 ^) 2 



cos 6i 



(4.22) 



The probability of detecting a muon decay at B is then: 

lDecm n Y ' n 



P{e + W\L,6 1 ) = 



E n 
x exp 

x/l. 



< e + u u \T R \fi + > 



cos #i 
sin 29 cos 



< u fi + \T R \n + > 

2 



2m 2 Am 2 



E 2 (ml 



ml) 2 ) _ 



m;E n 



- mlE u )L 



cos 6i 



(4.23) 
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where loec is defined in the same way as in Eqn(4.15). 



5 Discussion 

The quantum mechanics of neutrino oscillations has been surveyed in recent review 
articles [22, 23, 24], where further extensive lists of references may be found. 

In this Section, the essential differences between the calculations presented in the 
present paper and all previous treatments in the literature of the quantum mechanics of 
neutrino oscillations, as cited in the above review articles, will be summarised. A critical 
review of the existing literature will then be given. 

Hitherto, it has been assumed that the neutrino source produces a 'lepton flavour 
eigenstate' that is a superposition of mass eigenstates, at some fixed time. In this paper 
it is, instead, assumed following Shrock [25, 26] that the neutrino mass eigenstates are 
produced incoherently in different physical processes. This follows from the structure of 
the leptonic charged current in the Electroweak Standard Model: 

j^ cc ye P t = £^ a7/i (i _ l5 )u m ^ (5.1) 

Only the wavefunctions of the physical neutrino mass eigenstates appear in this cur- 
rent, and hence in the initial or final states of any physical process. A consequence is that 
the neutrino mass eigenstates can be produced at different times in the path amplitudes 
corresponding to different mass eigenstates. It has recently been shown that experimental 
measurements of the decay width ratio; r(7r — > eis)/T(ir — > fiu) and of the MNS matrix 
elements are inconsistent with the production of a coherent 'lepton flavour eigenstate' in 
pion decay [27] and that the the 'equal time' or 'equal velocity' hypothesis resulting from 
this assumption underestimates, by a factor of two, the contribution of neutrino propa- 
gation to the oscillation phase [28]. As demonstrated above, allowing for the possibility 
of different production times of the neutrinos results in an important, decay process de- 
pendent, contribution to the oscillation phase from the propagator of the source particle. 
The non-diagonal elements of U a i in Eqn(5.1) describe violation of lepton flavour (or gen- 
eration number) by the weak charged current. For massless neutrinos, the MNS matrix 
becomes diagonal; lepton flavour is conserved within each generation, and the familiar 
'lepton flavour eigenstates' are give by the replacements: v 1 — > u e , v 2 — > v z — > v T . 
Only in this case are the lepton flavour eigenstates physical, being all mass eigenstates of 
vanishing mass. 

The standard derivation of the neutrino oscillation phase will now be considered, 
following the treatment of Ref. [3], but using the notation of the present paper. The 
calculation is performed assuming an initial 'lepton lavour eigenstate of the neutrino' 
that is a superposition of the mass eigenstates \v\ > and \v 2 > : 

\vn >= U^i\vi,p > +U^ 2 \v 2) p >= - sin 9\u 1: p > + cos6\u 2 ,p > (5.2) 

where \vi,p > are mass eigenstates of fixed momentum p. This flavour eigenstate is 
assumed to evolve with laboratory time, t, according to fixed energy solutions of the 
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non-relativistic Schrodinger equation into the mixed flavour state \a,t >: 

\a,t >= - sm9e~ lElt \v 1: p > + cos 9e~ lE2t \is 2 ,p > (5.3) 

where Ei, E 2 are the laboratory energies of the neutrino mass eigenstates. The amplitude 
for transition into the 'electron flavour eigenstate': 

\u e ,p >= U e i\vi,p > +U e2 \u 2 ,p >= cos9\ui,p > + sm9\u 2 ,p > (5.4) 

at time t is then, using Eqns(5.3), (5.4): 

< i/ e , p\a, t >= sin 9 cos 9 (-e~ iElt + e~ iE ' 2t ) (5.5) 

Because it is assumed that the neutrinos have the same momentum but different energies : 

2 

I 777 

E t = vV + mf = P+^ + 0(mt) (5.6) 

and using (5.5) and (5.6), the probability of the flavour state v e at time t is found to be: 

P(u e , t) = | < v e , p\a, t>\ 2 = 2 cos 2 9 sin 2 9 ^1 + cos 

Finally, since the velocity difference of the neutrino mass eigenstates is 0(Am 2 ), then, to 
the same order in the oscillation phase, the replacement t — > L can be made in Eqn(5.7) 
to yield the standard oscillation phase of Eqn(l.l). 

The following comments may be made on this derivation: 

(i) The time evolution of the neutrino mass eigenstates in Eqn(5.3) according to the 
Schrodinger equation yields a non-Lorentz-invariant phase ~ Et, to be compared 
with the Lorentz-invariant phase ~ m 2 t/E given in Eqn(2.14) above. Although the 
two expressions agree in the non-relativistic limit E ~ m it is clearly inappropriate 
to use this limit for the description of neutrino oscillation experiments. It may be 
noted that the Lorentz-invariant phase is robust relative to different kinematical 
approximations. The same result is obtained to order m 2 for the phase of spatial 
oscillations independent of whether the neutrinos are assumed to have equal mo- 
menta or energies. This is not true in the non-relativistic limit. Assuming equal 
momenta gives the standard result of Eqn(l.l), whereas the equal energy hypothe- 
sis results in a vanishing oscillation phase. A contrast may be noted here with the 
standard treatment of neutral kaon oscillations, which follows closely the derivation 
in Eqns(5.2) to (5.7) above, except that the particle phases are assumed to evolve 
with time according, to the Lorentz invariant expression, exp[— imr], where m is 
the particle mass and r is its proper time, in agreement with Eqn(2.11). 

(ii) As pointed out in Ref. [6], the different neutrino mass eigenstates do not have 
equal momenta as assumed in Eqns(5.2) and (5.6). The approximation of assuming 
equal momenta might be justified if the fractional change in the momentum of 
the neutrino due to a non-vanishing mass were much less than that of the energy. 
However, in the case of pion decay as is readily shown from Eqns(2.4) and (2.6) 
above, the ratio of the fractional shift in momentum to that in energy is actually 
(m 2 + m 2 )/(m 2 — m 2 ) = 3.67; so, in fact, the opposite is the case. 



2p 



(5.7) 
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(iii) The derivation of Eqn(5.7) is carried out in the abstract Hilbert space of the neu- 
trino mass eigenstates or 'lepton flavour eigenstates' without any reference to the 
production or detection processes necessary for the complete description of an ex- 
periment in which 'neutrino oscillations' may be observed. In this calculation the 
'mass' and 'flavour' bases are treated as physically equivalent. However in Standard 
Model amplitudes only states of the mass basis appear. Also it has been pointed out 
that 'flavour momentum eigenstates' cannot be defined in a theoretically consistent 
manner [29]. Their existence is, in any case, excluded by experiment for the case of 
pion decay [27]. 

(iv) What are the physical meanings of t, p in Eqn(5.3)? In this equation it is assumed 
that the neutrino mass eigenstates are both produced, and both detected, at the 
same times. Thus both have the same time-of-flight t. The momentum p cannot be 
the same for both eigenstates, as assumed in Eqn(5.6), if both energy and momentum 
are conserved in the decay process. For any given value of the laboratory time t 
the different neutrino mass eigenstates must be at different space-time positions 
because they have different velocities 13 , if it is assumed that both mass eigenstates 
are produced at the same time. It then follows that the different mass eigenstates 
cannot be probed, at the same space-time point, by a neutrino interaction, whereas 
the latter must occur at a definite space-time point in every detection event. In fact, 
there is an inconsistent treatment of the velocity of the neutrinos. Equal production 
times imply equal space-time velocities, whereas it is assumed thast 'kinematical 
velocities' defined as Pi/Ei are different for the different mass eigenstates. 

(v) The historical development of the calculation of the neutrino oscillation phase is of 
some interest. The first published prediction [30] actually obtained a phase a factor 
two larger than Eqn(l.l) i.e. in agreement with the contribution from neutrino 
propagation found in the present paper. This prediction was later used, for example, 
in Ref. [31] . The derivation sketched above, leading to the standard result of Eqn(l. 1) 
was later given in Ref. [32]. A subsequent paper [33] by the authors of Ref. [31], 
published shortly afterwards, cited both Ref. [30] and Ref. [32], but used now the 
prediction of the latter paper. No comment was made on the factor of two difference 
in the two calculations. In a later review article, [34], by the authors of Ref. [31] 
a calculation similiar to that of Ref. [32] was presented in detail. Subsequently, 
all neutrino oscillation experiments have been analysed on the assumption of the 
standard oscillation phase of Eqn(l.l). 

It may be thought that the kinematical and geometrical inconsistencies mentioned in 
points (ii) and (iv) above result from a too classical approach to the problem. After all, 
what does it mean, in quantum mechanics, to talk about the 'position' and 'velocity' of 
a particle, in view of the Heisenberg Uncertainty Relations [35]? This point will become 
clear later in the present discussion, but first, following the original suggestion of Ref. [36], 
and, as done in almost all subsequent work on the quantum mechanics of neutrino oscilla- 
tions, the 'wave packet' description of the neutrino mass eigenstates will be considered. In 
this approach, both the 'source' and also possibly the 'detector' in the neutrino oscillation 
experiment are described by coherent spatial wave packets. Here the 'source' wavepacket 

13 This is true not only in the case of energy-momentum conservation, but also if it is assumed, as in the 
derivation of the standard formula, that the neutrinos have the same momentum but different energies. 
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treatment in the covariant approach of Ref. [7] will be briefly sketched. After discussing 
the results obtained, and comparing them with those of the present paper, the general 
consistency of the wave packet approach with the fundamental quantum mechanical for- 
mula (1.2) will be examined. A further discussion of wave packets as applied to neutrino 
oscillations can be found in Ref. 



The basic idea of the wave packet approach of Ref. [7] is to replace the neutrino prop- 
agator (2.11) in the path amplitude by a four-dimensional convolution of the propagator 
with a 'source wave packet' which, presumably, describes the space-time position of the 
decaying pion. For mathematical convenience, this wave packet is taken to have a Gaus- 
sian form with spatial and temporal widths a x and a t respectively. Thus, the neutrino 
propagator D is replaced by D where: 



D(x f - Xi,mj) = J d 4 xD(x f - x, ,mj)ip in (x - x { ,j) 



(5.8) 



where X{ and Xf are 4- vectors that specify the neutrino production and detection positions, 
respectively, and: 



ip in (x - Xi, j) = iVoexp 



-ipj ■ (x - Xi) - 



(x-Xi) 2 (t-uy 



4a? 



(5.9) 



The integral in (5.8) was performed by the stationary phase method, yielding the result 
(up to multiplicative and particle flux factors), and here assuming, for simplicity, one 
dimensional spatial geometry: 



D(Ax, At, rrij) = exp 



7X1 % 



-i\-fAt-P 3 (Ax-v 3 At)\- 



(Ax - VjAt) 2 



(5.10) 



where Ax = Xf — Xi and At — tf — U. For the case of — > v e oscillations', following n + 
decay at rest, the probability to observe flavour v e at time t and distance x is given by: 



P{v e , Ax, At) = cos 2 9 sin 2 9 -D(Ax, At, mi) + D(Ax, At, m 2 ) 



(5.11) 



Performing the integral over At and making the ultra-relativistic approximation v±, v-i ~ 1 
yields finally, with Ax = L: 



P(Ve,L) 



2 cos 2 9 sin 2 6(1 + 
exp 



Am\al + a 2 ) 
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(5.12) 



It can be seen that the oscillation phase in Eqn(5.12) is the same as standard one of 
Eqn(l.l). This is a consequence of the 'equal production time' hypothesis implicit in 
Eqn(5.10), where At does not depend on the mass eigenstate label j. The exponential 
damping factors in Eqn(5.12) are the same as those originally found in Ref. [37] for spatial 
wave packets (a t = 0). Considering now only spatial wave packets and using the property 
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a x a v = 1/2 derived from the Fourier transform of a Gaussian, the two terms in the 
exponential damping factor may be written as: 

Am 4 



F p = exp 

and 



(5.13) 



F x = exp 



Am 4 (m 2 +m 2 ) 2 L 
4« - rr^faf 



X J 



(5.14) 



The spatial damping factor F x is usually interpreted in terms of a 'coherence length' 

If 

2(m 2 -m 2 )V 



L » ~a"Z fir £ ( 5 - 15 ) 

then F x <C 1 and the neutrino oscillation term is strongly suppressed. Eqn(5.15) expresses 
the condition that oscillations are only observed provided that the wave packets overlap. 
Since Av ~ Am 2 the separation of the wave packets is ~ Av L ~ Am 2 L, so that Eqn(5.15) 
is equivalent to AvL ^> u x (no wave packet overlap). 

The damping factor F p is typically interpreted [39] in terms of the 'Heisenberg Un- 
certainty Principle'. This factor is small, unless the difference in mass of the eigenstates 
is much less than 2^ ^/m^cTp, so it is argued that only for wide momentum wave pack- 
ets can neutrino oscillations be observed, whereas in the contrary case, when the mass 
eigenstates are distinguishable, the interference effect vanishes. In the case of pion decay 
the difference in momentum of the two interfering mass eigenstates comes only from the 
Si term in Eqn(2.4), as 5 n , being a property of the common initial state, is the same for 
both eigenstates. The neutrino momentum smearing in pion decay is then estimated from 
Eqn(2.4) as: 



2Pn8jm,i r„m„ 

__ u 1 fj H V - ill ' in- 

m 2 — 2m n 



For Am 2 = l(eV) 2 the value of F p is found to be exp(— 500) ~ 7.1 x 10~ 218 giving complete 
suppression of neutrino oscillations. This prediction is in clear contradiction with the tiny 
damping corrections found in the path amplitude analysis in of Section 2 above. 

The preceding discussion of the derivation of the standard formula for the oscillation 
phase (1.1) in terms of 'flavour eigenstates' revealed contradictions and inconsistencies if 
the neutrinos are assumed to follow classical space-time trajectories. The 'source wave 
packet' treatment gives the standard result for the oscillation phase and predicts that the 
interference term will be more or less damped depending on the widths in space-time and 
momentum space of the wave packets. So do wave packets actually play a role in the 
correct quantum mechanical description of neutrino oscillations, as suggested in Ref. [36]? 
Are the packets actually constrained by the Heisenberg Uncertainty Relations? How do 
the properties of the wave packet effect the possibility to observe neutrino oscillations? 
The answers to all these questions are contained in the results of the calculations pre- 
sented in Section 2 above. They are now reviewed, with special emphasis on the basic 
assumptions made and the physical interpretation of the equations. 

Referring again to Fig.l, in a) a single pion comes to rest in the stopping target T. 
The time of its passage is recorded by the counter Ca, which thus defines the initial 
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Table 1: Different contributions to neutrino momentum smearing in pion decay at rest 
(see text). 

state as a n + at rest at time to- This pion, being an unstable particle, has a mass W 7 ^ 
that is, in general, different from its most likely value which is the pole mass m n . What 
are shown in Fig. lb) and Fig.lc) are two different classical histories of this very same 
pion. In b) it decays into the neutrino mass eigenstate v\ at time t±, and in c) it decays 
into the neutrino mass eigenstate z/ 2 at time t 2 . Because these are independent classical 
histories, the physical masses, and W^ 2 ), of the muons recoiling against the mass 

eigenstates v 1 and u 2 , respectively, are, in general, not equal. Taking into account, now, 
exact energy-momentum conservation in the decay processes (appropriate because of the 
covariant formulation used throughout) the eigenstates v\ and v 2 will have momenta which 
depend on W n and W^m and W n and W^) respectively. These momenta are calculated in 
Eqn(2.4). The distributions of W n and are determined by Breit-Wigner amplitudes 
(that are the Fourier transforms of the corresponding exponential decay laws) in terms 
of the decay widths T n and respectively. In accordance with Eqn(1.2), only the Breit- 
Wigner amplitudes corresponding to the physical muon masses are (coherently) integrated 
over at the amplitude level. The integral over (a property of the initial state) is 
performed (incoherently) at the level of the transition probability. Because of the long 
lifetimes of the n and //, the corresponding momentum wave packets are very narrow, 
so that all corrections resulting from integration over the resulting momentum spectra 
are found to give vanishingly small corrections. Indeed, as is evident from the discussion 
of Eqn(5.10) above, the width of the momentum wave packet is much smaller than the 
difference in the momenta of the eigenstates expected for experimentally interesting values 
of the neutrino mass difference (say, Am 2 = l(eV) 2 ). However, contrary to the prediction 
of Eqn(5.10), this does not at all prevent the observation of neutrino oscillations. This 
is because the oscillations result from interference between amplitudes corresponding to 
different propagation times, not different momenta, of the neutrino mass eigenstates. 
Table 1 shows contributions to the fractional smearing of the neutrino momentum from 
different sources: 

(1) neutrino mass difference Am 2 = l(eV) 2 (Eqn(2.4)) 

(2) coherent effect due to the physical muon mass (Si = T^/2 in Eqn(2.4)) 

(3) incoherent effect due to the physical pion mass (5 n = in Eqn(2.4)) 

(4) incoherent Doppler effect assuming p n = 2.6 x lCT 3 MeV 

The pion mass effect is of the same order of magnitude as the neutrino mass shift. The 
muon mass effect is two orders of magnitude smaller, while the Doppler effect at room 
temperature gives a shift eleven orders of magnitude larger than a (leV) 2 neutrino mass 
difference squared. 

According to the usual interpretation of the Heisenberg Uncertainty Principle, the 
neutrinos from pion decay, which, as has been shown above, correspond to very narrow 
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momentum wave packets, would be expected to have a very large spatial uncertainty. 
Indeed, interpreting the width of the coherent momentum wave packet generated by the 
spread in according to the the momentum-space Uncertainty Relation ApAx = 1 
gives Ax = 1.27 km. Does this accurately represent the knowledge of the position of a 
decay neutrino obtainable in the experiment shown in Fig.l? Without any experimental 
difficulty, the decay time of the pion can be measured with a precision of 10~ 10 sec, by 
detecting the decay muon. Thus, at any later time, the distance 14 of the neutrino from 
the decay point is known with a precision of c x l(U 10 cm = 3cm. This is a factor 4xl0 4 
more precise than the 'uncertainty' given by the Heisenberg relation. It is clear that the 
experimental knowledge obtainable on the position of the neutrino is essentially classical, 
in agreement with the theoretical description in terms of classical particle trajectories 
in space-time. In this case the momentum-position Uncertainty Relation evidently does 
not reflect the possible experimental knowledge of the position and momentum of the 
neutrino. This is because it does not take into account the prior knowledge that the mass 
of the neutrino is much less than that of the pion or muon, so that its velocity is, with 
negligible uncertainty, c. In spite of this, a Heisenberg Uncertainty Relation is indeed 
respected in the pion decay process. The Breit-Wigner amplitude that determines the 
coherent spread of neutrino momentum is just the Fourier transform of the exponential 
decay law of the muon. The width parameter of the Breit-Wigner amplitude and the 
muon mean lifetime do indeed respect the energy-time Uncertainty Relation T^r^ = I. 

It is then clear, from this careful analysis of neutrino oscillations following pion decay 
at rest, that, in contradiction to what has been almost universally assumed until now, 
the neutrinos are not described by a coherent spatial wave packet. There is a coherent 
momentum wave packet, but it is only a kinematical consequence of a Breit Wigner 
amplitude. If for mathematical convenience, the momentum wave packet is represented 
by a Gaussian, a conjugate (and spurious) Gaussian spatial wave packet will be generated 
by Fourier transformation. Indeed, in the majority of wave packet treatments that have 
appeared in the literature, Gaussian momentum and spatial wave packets related by a 
Fourier transform with widths satisfying the 'Uncertainty Relation' o v a x = 1/2 have 
been introduced. In fact, it is evident, by inspection, that the space-time wavepacket of 
Eqn(5.10) does not correctly reflect the known space-time structure of the sequence of 
events corresponding to pion decay. Since Ax = L, the fixed source-detector distance, 
the integral of At over the range from — oo to +oo assumed in order to derive Eqn(5.12) 
implies that the neutrino velocities v v = Ax/ At vary also in the range: — oo < v v < +oo. 
This unphysical range of integration results in an average neutrino momentum that is less 
than the kinematical velocity, Pi/Ei, of either mass eigenstate, due to contributions to the 
integral from negative values of At [27]. Also the production time dependence is known 
to be exponential, not Gaussian. There are indeed 'Heisenberg Uncertainties' in the 
production times of the neutrinos, due to the finite source lifetime, but once the neutrinos 
are produced their motion in space-time is well approximated, in each alternative history, 
by that of a free classical particle. The above discussion shows clearly that the ad hoc 
Gaussian wave packet introduced in Eqns(5.8) and (5.9) does not correspond to the actual 
sequence of the space-time events that constitute realistic neutrino oscillation experiments 

14 Here, 'distance' is defined without specifying the flight direction of the neutrino. Precise, simultane- 
ous, measurement of the flight direction of the recoil muon also determines, with a similar precision, the 
neutrino direction. Evidently the thickness of the stopping target may be chosen sufficiently small that 
its contribution to the neutrino position uncertainty is negligible. 
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The limitation on the detection distance L, for observation of neutrino oscillations, 
given by the damping factor F^fT^) of Eqn(2.30) is easily understood in terms of the 
classical particle trajectories shown in Fig.l. For a given velocity difference, the time 
difference t 2 — t 1 becomes very large when both neutrinos, in the alternative classical 
histories, are required to arrive simultaneously at a far distant detector. Because of the 
finite pion lifetime, however, the amplitude for pion decay at time t 2 is smaller than that 
at t\ by the factor exp[— (t 2 — ti)T 7r /2]. Integrating over all decay times results in the 
exponential damping factor F^rV) of Eqn(2.30). It is clear that, contrary to the damping 
factor F x of Eqn(5.14), the physical origin of the L dependent damping factor is quite 
unrelated to 'wavepacket overlap'. 

For a given value of t 2 — t\ the coherent neutrino momentum spread originating in the 
Breit Wigner amplitude for produces a corresponding velocity smearing that reduces 
the number of possible classical trajectories arriving at the detection event. This effect is 
taken into account in the integral shown in Eqn(2.22). The effect is shown, in Section 2 
above, to be much smaller than the (already tiny) pion lifetime damping described above, 
and it is neglected in Eqn(2.26). As described by Eqn(2.25), the (incoherent) integration 
over the Breit Wigner function containing W n gives an additional damping correction to 
the interference term that is also very tiny compared to that due to the pion lifetime. 

A final remark is now made on the physical interpretation of the damping factors (5.13) 
and (5.14) that have often been derived and discussed in the literature. F x is replaced, 
in the path amplitude calculations, by F U (T W ) and F p by the factor resulting from the 
coherent integration over the physical muon mass W^. The reason for the huge suppression 
factor predicted by Eqn(5.13), and the tiny one found in the path amplitude calculation, 
is that, in deriving (5.13) and (5.14), it is assumed that the neutrino eigenstates are both 
produced and detected at equal times. This will only be possible if both the hypothetical 
'wave packet overlap' is appreciable (Eqn(5.14)) and also the momentum smearing is 
sufficiently large that the time-of-flight differences due to the different neutrino masses 
are washed out (Eqn(5.13). In the path amplitude calculation interference and hence 
oscillations are made possible by different decay times of the source pion, and the damping 
factors analagous to F x and F p turn out to give vanishingly small corrections to the 
oscillation term. 

It may be remarked that the physical interpretation of 'neutrino oscillations' provided 
by the path amplitude description is different from the conventional one in terms of 
'flavour eigenstates'. In the latter the amplitudes of different flavours in the neutrino are 
supposed to vary harmonically as a function of time. This may be done, for example, by 
changing the basis states, in Eqn(5.3) above, from the mass to the flavour basis by using 
the inverses of Eqns(5.2) and (5.4). In the amplitudes for the different physical processes 
in the path amplitude treatment there is, instead, no variation of the 'lepton flavour' in 
the propagating neutrinos. If the mass eigenstates are represented as superpositions of 

15 It follows from this that discussions of quantum mechanical coherence in neutrino oscillations based 
on the properties of Gaussian wavepackets, as in Ref. [40], do not address the actual physical basis of 
neutrino oscillation experiments. A similar remark applies to the extensive study of Ref. [41], as well as 
the recent Ref. [42] in which it is claimed that spatial wave packets are necessary for a correct description 
of neutrino oscillations. 
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'lepton flavour eigenstates' there is evidently no temporal variation of the lepton flavour 
composition within each interfering amplitude. Indeed, the neutrinos in general occupy 
different spatial positions at any given time, making it impossible to project out a 'flavour 
eigenstate' at any time by using the inverses of Eqns(5.2) and (5.4). Only in the detection 
process itself where the different neutrino histories occupy the same space time point are 
the 'lepton flavour eigenstates' projected out, and the interference effect occurs that is 
described as 'neutrino oscillations'. In the case of the observation of decay products of 
the recoil muons no such projection on to a 'lepton flavour eigenstate' takes place, but 
exactly similar interference effects are predicted to occur. As previously emphasised [8], 
the 'flavour oscillations' of neutrinos, neutral kaons and b-mesons are just special examples 
of the universal phenomenon of quantum mechanical superposition described by Eqn(1.2), 
that also describes all the interference effects of physical optics. 

A more detailed critical review will now be made of previous treatments of the quantum 
mechanics of both neutrino and muon oscillations in the literature. 

By far the most widespread difference from the path amplitude treatment of the 
present paper is the non-respect of the basic quantum mechanical formula, (1.2), by the 
introduction of wave packets to describe 'source' and/or 'detector' particles[7,8, 10, 11,30- 
46]. Since in any practical neutrino oscillation experiment a single initial or final quantum 
state, as specified by Eqn(1.2), is not defined, but rather sets of initial and final states 
I — Hik an d F = J2 m fm determined by experimental conditions, Eqn(1.2) may be 
generalised to: 
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to be contrasted with the formula used in the references cited above 



(5.16) 
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Here, ipi and ipf are 'source' and 'detector' wave packets respectively. In Eqn(5.17) the 
initial states i\ and final states f m of Eqn(5.16), that correspond to different spatial posi- 
tions, and also, possibly, different kinematic properties, of the source particle or detection 
event, are convoluted with ad hoc spatial and/or temporal 'wave packets', that are, in 
the opinion of the present writer, for the reasons given above, devoid of any physical 
significance. 

A possible reason for the widespread use of Eqn(5.17) instead of (5.16) may be under- 
stood following a remark of the author of Ref. [46] concerning a 'paradox' of the complete 
quantum field theory calculation that takes into account, by a single invariant amplitude, 
production, propagation and detection of the neutrinos. It was noticed that, if the ampli- 
tude for the complete chain of processes is considered to correspond to one big Feynman 
diagram, then integration over the space-time coordinates of the initial and final states 
particles will reduce the exponential factors, containing the essential information on the 
interference phase, to energy-momentum conserving delta functions, and so no oscilla- 
tions will be possible. A related remark was made by the authors of Ref. [11] who stated 
that, as they were assuming exact energy-momentum conservation, the integration over 
the space-time coordinates could be omitted. They still, however, (quite inconsistently, in 
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view of the previous remark) retained the exponential factors containing the interference 
phase information. These considerations indicate a general confusion between momen- 
tum space Feynman diagram calculations, where it is indeed legitimate to integrate, at 
the amplitude level, over the unobserved space-time positions of the initial and final state 
particles, and the case of neutrino oscillations, where the amplitude is defined with initial 
and final states corresponding to space-time positions. In the latter case, it is the unob- 
served momenta of the propagating particles that should be integrated over, as is done in 
Eqns(2.1) and (2.35) above, and not the space-time positions of the 'source' or 'detector' 
particles, as in Eqn(5.17). 

It is clear that exact energy-momentum conservation plays a crucial role in the path 
amplitude calculation. This is valid only in a fully covariant theory. Still, several authors, 
in spite of the ultra-relativistic nature of neutrinos, used a non-relativistic theory to 
describe the production, propagation and detection of neutrinos [44, 47, 51, 52]. As is 
well known, in such 'Old Fashioned Perturbation Theory' [54] energy is not conserved 
at the level of propagators and so no precise analysis of the kinematics and the space- 
time configurations of the production and detection events, essential in the covariant path 
amplitude analysis, is possible. 

Even when, in some cases, the complete production, propagation and detection process 
of the neutrinos were described [41, 43, 44, 46, 47], equal neutrino energies [44, 46, 47], 
equal neutrino momenta [43] or either [41] were assumed, in contradiction with energy- 
momentum conservation and a consistent space-time description of the production and 
detection events. As follows directly from Eqn(5.6) (or the similar formula, for the neu- 
trino momentum, obtained by assuming equal neutrino energies), the standard formula 
(1.1) for the oscillation phase was obtained in all the above cited references. As shown 
in Reference [28] this is a consequence of the universal equal production time (or equal 
velocity) hypothesis. In fact, the assumptions of equal momenta, equal energies or exact 
energy-momentum conservation give only negligible, 0(m 4 ), changes in the oscillation 
phase [28]. 

An interesting discussion of the interplay between different kinematical assumptions 
(not respecting energy-momentum conservation) and the space-time description of the 
production and detection events was provided in Ref. [5]. This treatment was based 
on the Lorentz-invariant propagator phase of Eqn(2.11). By assuming either equal mo- 
mentum or equal energy for the propagating neutrinos, but allowing different times of 
propagation for the two mass eigenstates, values of <$>\ 2 agreeing with Eqns(2.14) and 
(2.24) above were found, i.e. differing by a factor two from the standard formula. Al- 
teratively, assuming equal velocities, (and hence equal propagation times) the standard 
result (1.1) was obtained. In this latter case, however, the masses, momenta and energies 
of the neutrinos must be related, up to corrections of 0(m|), according to: 

- = § = f (5-18) 
m 2 P 2 E 2 K ' 

Since the ratio of the neutrino masses may take, in general, any value, so must then the 
ratio of their momenta. For the case of neutrino production from pion decay at rest with 
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mi, m2 "C m^,m n the relation (5.18) is clearly incompatible with Eqn(2.4) which gives: 

^ - 1 - Am2K + m ' } + Ofm 4 m 4 ) f5 19) 

P 2 ~ L K-m^ + <Am 1 ,m 2 j (5-19) 

On the other hand, Eqn(5.19) is clearly compatible (up to corrections of 0(Am 2 )) with 
the 'equal momentum' hypothesis. There is similar compatiblity with the 'equal energy' 
hypothesis. Even so, the authors of Ref. [5] recommended the use of the equal velocity 
hypothesis. With the hindsight provided by the path amplitude analysis, in which the 
two neutrino mass eigenstates do indeed have different propagation times, it can be seen 
that the kinematically consistent 'equal momentum' and 'equal energy' choices are good 
approximations and the neutrino oscillation phase, resulting from the propagation of the 
neutrinos alone, is indeed a factor two larger than the prediction of the standard formula. 
A more detailled discussion of the effects of the 'equal momentum' 'equal energy' and 
'equal velocity' hypotheses may be found in Reference [28] 

The author of Ref. [46] included the propagator of the decaying pion in the complete 
production-propagation-detection amplitude; compare Eqn(8) of Ref. [46] with Eqn(2.1) 
above. However, no detailed space-time analysis of production and detection events was 
performed. Square spatial wave packets for the 'source' and 'detector' were convoluted 
at amplitude level as in Eqn(5.17). As the neutrinos were assumed to have a common 
production time, no contribution to the oscillation phase from the source particle was 
possible, and so the standard result for the oscillation phase was obtained. 

Although, following Ref. [44], most recent studies of the quantum mechanics of neu- 
trino oscillations have considered the complete production-propagation-detection process, 
some authors still use, in spite of the criticisms of Ref. [29], the 'flavour eigenstate' de- 
scription [4, 55, 56]. In the last two of these references a 'quantum field theory' approach 
ia adopted, leading to an 'exact' oscillation formula [56] that does not make use of the 
usual ultra-relativistic approximation. This work included neither exact kinematics, nor 
an analysis of the space-time structure of production and detection events. In Ref. [4] 
the equal energy hypothesis was used and in Ref. [55, 56] the equal momentum hypoth- 
esis. In all three cases the standard result was found for the oscillation phase in the 
ultra-relativistic limit, as a consequence of the assumption of equal production times. 

Correlated production and detection of neutrinos and muons produced in pion decay 
were considered in Ref. [8]. The introduction to this paper contains a valuable discussion 
of the universality of the 'particle oscillation' phenomenon. It is pointed out that this is 
a consequence of the general principle of amplitude superposition in quantum mechanics, 
and so is not a special property of the K° — K°, B° — B° and neutrino systems which are 
usually discussed in this context. This paper used a covariant formalism that employed the 
'energy representation' of the space-time propagator. In the introduction, the important 
difference between Eqns(5.16) and (5.17) was also touched upon: 

' The reader will agree that one should not integrate over space if one is interested in 
spatial interference (or oscillation).' 

Even so, in the amplitude for the correlated detection of the muon and neutrino 
(Eqn(2.10) of Ref. [8]) not only are the space-time positions of the production points of 
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the neutrino and muon integrated over, but they are assumed to be at different space- 
time points. The propagator of the decaying pion is not included in the amplitude, and 
although exact energy-momentum conservation is imposed, no space-time analysis of the 
production and decay points is performed. Correlated spatial oscillations of neutrinos and 
muons are predicted, though with interference phases different from the results of both the 
present paper and the standard formula. Pion and muon lifetime effects were mentioned 
in Ref . [8] , but neither the role of the pion lifetime in enabling different propagation times 
for the neutrinos nor the momentum smearing, induced by the Fourier-transform-related 
Breit Wigner amplitudes, were discussed. 

The claim of Ref. [8]) that correlated neutrino-muon oscillations should be observable 
in pion decay was questioned in Ref. [11]. The authors of the latter paper attempted to 
draw conclusions on the possibility, or otherwise, of particle oscillations by using 'plane 
waves', i.e. energy-momentum eigenfunctions. As is well known, such wavefunctions are 
not square integrable, and so can yield no spatial information. The probability to find a 
particle described by such a wave function in any finite spatial volume is zero. Due to the 
omission of the (infinite) normalisation constants of the wavefunctions many of the equa- 
tions in Ref. [11] are, as previously pointed out [57], dimensionally incorrect. Momentum 
wavepackets for the decaying pion convoluted at amplitude level as in Eqn(5.17) were 
also discussed in Ref. [11]. Although exact energy-momentum conservation constraints 
were used, it was assumed, as in Ref. [8], that the muons and the different neutrino mass 
eigenstates are both produced and detected at common points (Eqn(35) of Ref. [11]). The 
latter assumption implies equal velocities, yielding the standard neutrino oscillation phase 
as well as the inconsistent kinematical relation (5.18). The authors of Ref. [11] concluded 
that: 

(a) correlated /i — v oscillations of the type discussed in Ref. [8] could be observed, 
though with different oscillation phases. 

(b) oscillations would not be observed if only the muon is detected 

(c) neutrino oscillations can be observed even if the muon is not detected. 

Conclusion (b) is a correct consequence of the (incorrect) assumption that the muons 
recoiling against the different neutrino mass eigenstates have the same velocity. As both 
muons have the same mass they will have equal proper time increments. So according to 
Eqn(2.12) the phase increments will also be equal and the interference term will vanish. 
The conclusion (c) is in agreement with the prediction of Eqn(3.22) of Ref. [8]. The path 
amplitude calculation of the present paper shows that conclusion (b) is no longer valid 
when the different possible times of propagation of the recoiling muons are taken into 
account. 

Observation of neutrino oscillations following pion decay, using a covariant formalism 
(Schwinger's parametric integral representation of the space-time propagator) was con- 
sidered in Ref. [48]. Exact energy- momentum conservation was imposed, and integration 
over the pion spatial position at amplitude level, as in Eqn(5.17), was done. The propaga- 
tor of the pion source was included in the amplitudes, but as the different mass eigenstates 
were produced and detected at the same space-time points, equal propagation velocities 
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were implicitly assumed, so that just as in Refs[5, 11], where the same assumption was 
made, the standard neutrino oscillation phase was obtained. In the conclusion of this 
paper the almost classical nature of the space-time trajectories followed by the neutrinos 
was stressed, although this was not taken to its logical conclusion in the previous discus- 
sion, e.g. the kinematical inconsistency of the equal velocity hypothesis that requires the 
evidently impossible condition (5.18) to be satisfied. 

In a recent paper [58] , the standard neutrino oscillation formula with oscillation phase 
given by Eqn(l.l) was compared with a neutrino decoherence model. In order to take into 
account incertainties in the position of the source and the neutrino energy, an average was 
made over the quantity L/4E U , assuming that it is distributed according to a Gaussian 
with mean value i and width a. The average was performed in an incoherent manner. 
Thus the calculation is closely analagous to those for the effects of target or detector length 
or of thermal motion of the neutrino source, presented in the Appendix A of the present 
paper. Perhaps uniquely then, in the published literature, in Ref. [58] the effects of source 
position and motion are taken into account correctly, according to Eqn(5.16) instead of 
Eqn(5.17). However, the source of the neutrino energy uncertainty is not specified. In as 
far as it is generated from source motion the calculation is, in principle, correct. There 
is however also the (typically much smaller, see Table 1 above, for the case of pion decay 
at rest) coherent contribution originating from the variation in the physical masses of the 
unstable particles produced in association with the neutrino, as discussed in detail above. 
It was concluded in Ref. [58] that the Gaussian averaging procedure used gave equivalent 
results to the decoherence model for a suitable choice of parameters. 

It is clearly of great interest to apply the calculational method developed in the present 
paper to the case of neutral kaon and b-meson oscillations. Indeed the use of the invari- 
ant path amplitude formalism has previously been recommended [59] for experiments 
involving correlated pairs of neutral kaons. Here, just a few remarks will be made on the 
main differences to be expected from the case of neutrino or muon oscillations. A fur- 
ther discussion can be found in Ref. [28] and a more detailed treatment will be presented 
elsewhere [60]. 

In the case of neutrino and muon oscillations, the interference effect is possible as 
the different neutrino eigenstates can be produced at different times. This is because the 
decay lifetimes of all interesting sources (pions, muons, /3-decaying nuclei) are much longer 
than the time difference beween the paths corresponding to the interfering amplitudes. To 
see if a similar situation holds in the case of Ks — Kl oscillations, three specific examples 
will be considered with widely differing momenta of the neutral kaons: 

(I) - K S K L 

(II) n-p -> AK° at v/i = 2 GeV 
(III) n-p -> AK° at v/i = 10 GeV. 

These correspond to neutral kaon centre-of-mass momenta of 108 MeV, 750 MeV and 
5 GeV respectively. In each case the time difference (At K ) of production of K s and 
K L mesons, in order that they arrive at the same time at a point distant c^ K r s (where 
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7if is the usual relativistic parameter) from the source in the centre-of-mass frame is 
calculated. Exact relativistic kinematics is assumed and only leading terms in the mass 
difference Am^ = uil — ms are retained. Taking the value of Amx and the various 
particle masses from Ref. [61] the following results are found for Atx in the three cases: 
(I) 2.93xl0~ 24 sec, (II) 8.3xl0~ 25 sec and (III) 6.4x 10~ 26 sec. For comparison, for neutrino 
oscillations following pion decay at rest, with Am 2 = (leV) 2 and L = 30 m, Eqn(2.29) 
gives At u = 5.6 x lCT 23 sec. The result (I) may be compared with the mean life of the 
meson of 1.5xl0~ 22 sec [61]. Thus the <fi lifetime is a factor of about 27 larger than 
Atx indicating that Ks — Kl interference should be possible by a similar mechanism to 
neutrino oscillations following pion decays, i. e. without invoking velocity smearing of the 
neutral kaon mass eigenstates. In cases (I) and (II) the interference effects observed will 
depend on the 'characteristic time' of the non resonant (and hence incoherent) strong 
interaction process, a quantity that has, hitherto, not been susceptible to experimental 
investigation 16 . If this time is much less than, or comparable to, Atx, essentially equal 
velocities (and therefore appreciable velocity smearing) of the eigenstates will be necessary 
for interference to occur. Since Awlk and Ts are comparable in size, velocity smearing 
effects are expected to be, in any case, much larger than for neutrino oscillations following 
pion decay. These effects may be roughly estimated by using the Gaussian approximation 
(2.22) of the present paper. The main contribution to the velocity smearing is due to the 
variation of the physical mass of the K s rather than those of the K L or A. 

For the Bi — B 2 oscillation case, analagous to (I) above, Y(4S) — > BiB 2 , where pb = 
335 MeV, the value of Ats is found to be 1.4 x 10~ 22 sec, to be compared with t(T(AS)) = 
4.7 x 10~ 23 sec [61], which is a factor 3 smaller. Thus, velocity smearing effects are expected 
to play an important role in Bi — B 2 oscillations. This is possible, since the neutral b- 
meson decay width (4.3 x 1CT 10 MeV), and mass difference (3.1 x 10~ 10 MeV), have similar 
sizes. 

In closing, it is interesting to mention two types of atomic physics experiments where 
interference effects similar to the conjectured (and perhaps observed [17, 63, 64]) neutrino 
oscillations have aleady been clearly seen. 

The first is quantum beat spectroscopy [65]. This type of experiment, which has 
previously been discussed in connection with neutrino oscillations [57], corresponds closely 
to the gedanken experiment used by Heisenberg [15] to exemplify the the fundamental law 
of quantum mechanics, Eqn(1.2). The atoms of an atomic beam are excited by passage 
through a thin foil or a laser beam. The quantum phase of an atom with excitation energy 
E* evolves with time according to: exp(— iE*At) (see the discussion after Eqn(3.3) above). 
If decay photons from two nearby states with excitation energies E* a and Ejj are detected 
after a time interval At ( for example by placing a photon detector beside the beam at a 
variable distance d from the excitation foil) a cosine interference term with phase: 

(g - E* p )d 

(Pbeat — ~ 

^atom 

where v atom , is the average velocity of the atoms in the beam, is observed [65]. An atom 
in the beam, before excitation, corresponds to the neutrino source pion. The excitation 

16 A similar physical quantity has been considered in Ref. [62], where the possiblity of observable 
modifications to the exponential decay law and the Breit-Wigner line shape distribution is suggested. 
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(5.20) 



process corresponds to the decay of the pion. The propagation of the two different excited 
states, alternative histories of the initial atom, correspond to the alternative propagation 
of the two neutrino mass eigenstates. Finally the dexcitation of the atoms and the de- 
tection of a single photon corresponds to the neutrino detection process. The particular 
importance of this experiment for the path amplitude calculations presented in the present 
paper, is that it demonstrates, experimentally, the important contribution to the inter- 
ference phase of the space-time propagators of excited atoms, in direct analogy to the 
similar contributions of unstable pions, muons and nuclei discussed above. 

An even closer analogy to neutrino osillations following pion decay is provided by 
the recently observed process of photodetachement of an electron by laser excitation: 
the 'Photodetachment Microscope' [66]. A laser photon ejects the electron from, for 
example, an 16 0~ ion in a beam. The photodetached electron is emitted in an S-wave 
(isotropically) and with a fixed initial energy. It then moves in a constant, vertical, electric 
field that is perpendicular to the direction of the ion beam and almost parallel to the laser 
beam. An upward moving electron that is decelerated by the field eventually undergoes 
'reflection' before being accelerated towards a planar position-sensitive electron detector 
situated below the beam and perpendicular to the electric field direction (see Fig.l of 
Ref. [66]). In these circumstances, it can be shown [67] that, just two classical electron 
trajectories link the production point to any point in the kinematically allowed region of 
the detection plane. Typical parameters for 16 0~ are [68]: initial electron kinetic energy, 
102 /xeV: detector distance, 51.4 cm; average time-of-flight, 117 ns; difference in emission 
times for the electrons to arrive in spatial-temporal coincidence at the detector plane, 160 
ps. An interference pattern is generated by the phase difference between the amplitudes 
corresponding to the two allowed trajectories. The phase difference, derived by performing 
the Feynman path integral of the classical action along the classical trajectories [68], 
gives a very good description of the observed interference pattern. The extremely close 
analogy between this experiment and the neutrino oscillation experiments described in 
Sections 2 and 3 above is evident. Notice that the neutrinos, like the electrons in the 
photodetachement experiment, must be emitted at different times, in the alternative 
paths, for interference to be possible. This is the crucial point that was not understood 
in all previous treatments of the quantum mechanics of neutrino oscillations. 

Actually, Ref. [68] contains, in Section IV, a path amplitude calculation for electrons 
in free space that is geometrically identical to the discussion of pion decays in flight 
presented in Section 4 above (compare Fig.3 of the present paper with Fig. 3 of Ref. [68]) 
The conclusion of Ref. [68] is that, in this case, no interference effects are possible for 
electrons that are mononergetic in the source rest frame. As is shown in Section 4 above, 
if these electrons are replaced either by neutrinos of different masses from pion decay, or 
muons recoiling against such neutrinos, observable interference effects are indeed to be 
expected. 
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Appendix A 

Random thermal motion of the decaying pion in the target has two distinct physical 
effects on the phase of neutrino oscillations, 

012 (0) = p— + 2 (Al) 

(where the first and second terms in Eqn(Al) give the contributions of the neutrino and 
pion paths respectively): 

(1) The observed neutrino momentum, P u , is no longer equal to Po, due to the boost 
from the pion rest frame to the laboratory system. (Doppler effect or Lorentz boost) 

(2) The time increment of the pion path t D — t (see Eqn(2.16)) no longer corresponds 
to the pion proper time. (Relativistic time dilatation) 



Taking into account (1) and (2) gives, for the neutrino oscillation phase: 
where: 



Pu = 7tt-Po(1 + v n cos 0*) and ^ n 



Here, 0* is the angle between the neutrino momentum vector and the pion flight direction 
in the pion rest frame. Developing j w and v n in terms of the small quantity p^/m^, 
Eqn(A2) may be written as: 



012 \corr) = 0^(0) + 



Am 2 L 



1 - 



cos 9* 



Am 2 L 

2Pn 



1 - 



3m,,- 



2Pn 



(A3) 



Performing now the average of the interference term over the isotropic distribution in 

cos 6*: 



(cos 4>i2 {corr)) e * v 



—Re J exp [i(pi2 ( corT )] dcosOl 



^i?eexpii0^(O) + | — 



1 / Am 2 L 
6 [ P 





'(• 






1 - 






Po- 



Am 2 L 
2P 

2N 



1 - 



3m n 
~2Po~ 



{AA) 



In deriving Eqn(A4) the following approximate formula is used: 



e iac dc 



1 

2ia 



sin a a 

— — 

a o 



where 



a = 



Pk Am 2 L 
m w P 



m, 
~P< 



< 1 



(A5) 
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The average over the Maxwell-Boltzmann distribution (2.31) is readily performed by 'com- 
pleting the square' in the exponential, with the result: 



(cos 4>u (cor r))e t , P „ 



-Re exp 



'p n Am 2 L 
v 2m 7r P 



i _ Hh. 



+ i 



4 Km^J V P 



3m,r 
2P7 



(A6) 



= F"(Dpp) cos[0ET (0) + 0"(Dop)] 

leading to Eqns(2.32) and (2.33) for the Doppler damping factor F u (Dop) and phase shift 
4> u (Dop), respectively. 

The correction for the effect of thermal motion in the case of muon oscillations 
<fi v (Dop), respectively, is performed in a similar way. The oscillation phase: 

(A7) 



0^(0) = 



mlEgAm 2 L m^Am'L 



2m 7r P 3 



2P 3 



is modified by the Lorentz boost of the muon momentum and energy, and the relativistic 
time dilatation of the phase increment of the pion path, to: 



(j)t?(corr) 



mlE^Am 2 L m t ,Am 2 L 



(AS) 



where 



E ll = >y n ES(l + v n v8coa0*n) 
and Vq is given by Eqn(2.37). Developing, as above, in terms of p n /m n , gives: 



(f) 12 (cor r) 



^(0) + 



Pn VqW 2 Am 2 L 



'12 



^o 3 



n* i ( Pit 
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2 m 2 Am 2 L 
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^o 3 



2m T 



(A9) 



Performing the averages over #* and p w then leads to Eqns(2.53) and (2.54) for the damp- 
ing factor F^(Dop) and phase shift (f)^(Dop), respectively. 

The effect of the finite longitudinal dimensions of the target or detector is calculated 
by an appropriate weighting of the interference term according to the value of the distance 
X = Xf — Xi between the decay and detection points (see Fig.l). Writing the interference 
phase as i2 = /3X, and assuming a uniform distribution of decay points within the target 
of thickness £t- 



(cos0i 2 ) = — 



1 

It 
2 



sm ■ 



cos /3XdX 



cos/SL 



Ptr ' 2 

= F Targ cos(3L (A10) 

Substituting the value of (3 appropriate to neutrino oscillations yields Eqn(2.34). Since the 
value of (5 is the same for neutrino and muon oscillations, the same formula is also valid 
in the latter case. The same correction factor, with the replacement £t — > £d describes 
the effect of a finite detection region of length £ D : 



L — 



' D 



+ Xi < Xf < L + 



■D 



+ Xi 
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Appendix B 



The first step in the derivation of Eqn(4.17) relating Av(fi) to Am 2 is to calculate 
the angle 5*, in the centre-of-mass (CM) system of the decaying pion, corresponding to 
S in the laboratory (LAB) system (see Fig. 3). It is assumed, throughout, that the pion 
and muon are ultra-relativistic in the latter system, so that: v n ,Vi(fj) ~ 1. The Lorentz 
transformation relating the CM and LAB systems gives the relation: 

7 ff (l + v*(n) cos 9*) 

The starred quantities refer to the pion CM system. Making the substitutions: 9 2 = 0±+5, 
0* = 0*-\- S*, Eqns(Bl) may be solved to obtain, up to first order in 5, 5* and Am 2 : 

A,» = v*M - v*M = 7^(1 + v*M cosOtfS- vMicosdj + v*M)6* {m) 

sin 0^ 

where, (c.f. Eqn(2.37): 

m%. + m^ 

Using Eqn(2.36) Av*(fj) may be expressed in terms of the neutrino mass difference: 

. , 4m 2 m 2 Am 2 

Av » = m — V 2M 2 , 2 m 

Eliminating now Av*(/jl) between (B2) and (B4) gives a relation between 5, 5* and Am 2 : 

p = <y„(l + v*M cos61) 2 5 4m 2 m 2 Am 2 sin 6{ 

vt(n)(co&6t + vl(n)) (m 2 -m 2 ) 2 (m 2 +m 2 )(cos^ +v%(n)) 

In the LAB system, and in the UR limit, the difference of the velocities of the muons 
recoiling against the two neutrino mass eigenstates is: 

A„M = vM - vM = jgg - «g * |p 2 („) - EM] (BB) 

where is the muon energy in the LAB system for vanishing neutrino masses. Making 
the Lorentz transformation of the muon energy from the pion CM to the LAB frames, and 
using Eqns(2.4) and (2.36) to retain only terms linear in Am 2 and 5*, enables Eqn(B6) 
to be re-written as: 



AB( " ) = 4fe, 



Am 2 (cos^ ■ /•",(//)) 



Vo(fJ-) 



5*(m 7r — m M ) sin 9{ 



(B7) 



where E n is the energy of the pion beam. By combining the geometrical constraint 
equation for the muon velocities, (4.16) with (B5) and (B7) the angles 5 and 5* may be 
eliminated to yield the equation for LAB frame velocity difference: 

2m;(m; — m z ) B 
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where 



A = - u^cosfli) { (cos 6* +Vo(fx)) 2 + 



4mJ,mt sin 2 



9\ 



(m 



(59) 



x 



5 2 (m 2 + m 2 )(l + v*(fj.) cos 9l) 2 (cos 91 + ^(/i))(^(/x) - v n cos9 1 

m 2( m 2 _ m 2) 



4m 2 ,m 2 sin 2 0? 1 



2 + mlf 



To simplify (B8), the quantity (vi([/,)—v n cos 6 1 !) is now expressed in terms of kinematic 
quantities in the pion CM system. Within the UR approximation used, 

9 1 ,m 7T /E n ,m^/E^ < 1 

so that 

vM _ v , co S ^i(|-|| + + 0( , (^) * , «?) ( mi) 

Writing Eqn(Bl) to first order in 6> 1; and neglecting terms of 0(9 1 m 2 ): 

= m n v* sin 91 
1 E 7T (l + v*Mcos9* 1 ) 1 ; 

Using Eqn(B12), and expressing 5 M in terms of pion CM quantities, Eqn(Bll) may be 
written as: 

Expressing the RHS of (B13) in terms of E w and E^, using the relation: 

ml (2E U - E n ) - mlE n 

E^ml-ml) 1 J 

gives Eqn(4.20) of the text. 

On substituting (B13) into the RHS of (BIO), it can be seen that the factor in the 
large curly brackets is the same in (B9) and (BIO), and so cancels in the ratio A/B in 
Eqn(B8). It follows that: 

_ m 2 Am 2 / cosOl+vfo) \ 
™- E*(ml-m*) {l + v*Mcos9t) {B15) 

Finally, using (B3) and (B14) to express the factor in large brackets in Eqn(B15) in terms 
of E^ and E w , Eqn(4.17) of the text is obtained. 
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Figure 1: The space-time description of an experiment in which neutrinos produced in the 
processes n + — ► /i + (z^i, v 2 ) are detected at distance, L, via the processes (z/ 1; u 2 )n — > e~p. 
In a) a 7r + comes to rest in the stopping target T at time t . The pion, at rest at time t , 
constitutes the initial state for the path amplitudes. In b) and c) are shown two alternative 
classical histories for the 7r + ; in b),[ c)] the pion decays into the mass eigenstate \u x >, 
| 1 1/2 > ] at times ti, [ t 2 ]. If mi > m 2 , and for suitable values of t\ and t 2 (t 2 > h), 
the two classical histories may correspond to a common final state, shown in d) where 
the neutrino interaction (z/ 1; v 2 )n — > e~p occurs at time tp. As the initial and final states 
of the two classical histories are the same, the corresponding path amplitudes must be 
added coherently, as in Eqn(1.2), to calculate the probability of the whole process. 
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Figure 2: The space-time description of an experiment in which muons produced in 
the processes n + — > ii + {v ll u 2 ) are detected at distance, L, via decay processes /i + — > 
e + (ui, v 2 ){p\,v 2 ), denoted, conventionally, as — > e+z/eZ/^' As in Fig.l, b) and c) show 
alternative classical histories of the stopped tt + . If mi > m 2 the velocity of /xi is less 
than that of fi 2 , and provided that t 2 > tj., the muons may arrive at the same spatial 
point at the same time tjj in both classical histories. If the muons are detected at this 
space-time point in any way (not necessarily by the observation of muon decay as shown in 
c)) interference between the correponding path amplitudes occurs, according to Eqn(1.2), 
just as in the case of neutrino detection. 
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Figure 3: Two dimensional spatial geometry for the observation of neutrino or muon 
oscillations following pion decay in flight. Four possible classical histories of a pion, 
originally at the point A, are shown. In the first two, the pion decays either into the mass 
eigenstate \ui >, at point A or into \u 2 > at point C. If m 1 > m 2 , and for suitable values 
of the angles Q\ and 82, the neutrinos may arrive at the point B at the same time. If 
a neutrino detection event, such as (vi, 1/2)71 — > e~p, then occurs at B at this time, the 
paths AB and ACB will be indistinguishable so that the corresponding amplitudes must be 
superposed, as in Eqn(1.2), to calculate the probability of the overall decay-propagation- 
detection process. The third and fourth classical histories are similar, except that the 
neutrino mass eigenstates are replaced by the corresponding recoil muons. The muons in 
the different histories may arrive at point B, at the same time, leading to interference and 
'muon oscillations' if they are detected there. 



45 



